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SECTION  I 


INTRODUCTION 


PROGRAM  DESCRIPTION 


The  purpose  of  the  Structural  Life  Prediction  and  Analysis  Technology  research  and  devel- 
opment program  was  to  develop  and  verify  a systematic  method  for  predicting  fatigue  life 
exhaustion  of  military  gas  turbine  engine  disks  on  a mission  utilization  basis.  Major  emphasis 
was  given  to  the  evaluation  of  procedures  for  modeling  structural  response  to  mission  load- 
ing, the  development  of  a systematic  fatigue  life  exhaustion  specimen  test  program  and 
improved  fatigue  life  prediction  model,  and  verification  of  the  life  prediction  methodology 
with  full-scale  component  tests.  The  major  elements  of  a structural  life  prediction  method- 
ology are  shown  in  Figure  1 . 


The  program  was  divided  into  two  phases.  The  first  phase  concerned  fan  and  low-pressure 
compressor  disk  environments  where  dwell  (creep)  response  is  minimal;  the  second  phase 
concerned  high-pressure  compressor  and  turbine  environments  including  dwell  effects.  Full- 
scale  component  testing  for  each  phase  utilized  engine  disk  designs  which  contained  low 
cycle  fatigue  (LCF)  life  limited  bolt-holes  and  were  tested  to  representative,  complex  cycle, 
mission  histories.  The  first  phase  utilized  a prototypical  FI 00  fan  disk  geometry  and  titanium 
alloy  (Ti-6Al-2Sn-4Zr-6Mo)  tested  in  a servo-hydraulic  Ferris  Wheel  to  simulate  flight  loading 
conditions  at  the  bolt-holes.  The  second  pha^  utilized  a prototypical  advanced  transport 
high-pressure  turbine  disk  made  of  Waspaloy  , a nickel  base  superalloy. 

LCF  specimen  testing  was  performed  on  unnotched,  strain-controlled  (SC)  specimens, 
notched,  load-controlled  specimens  including  notched  round  bars  (NRB)  and  the  bolt-hole 
(BH)  subcomponent  specimen.  All  specimens  were  machined  from  the  same  heats  of  ma- 
terial used  to  obtain  the  full-scale  disks,  with  one  heat  for  each  alloy. 

The  nominal  (unnotched)  loading  conditions  for  both  phases  were  determined  from  mission 
history  analyses  performed  on  several  current  Air  Force  engine  models  using  a variety  of 
mission  conditions.  The  unnotched  loading  conditions  included  cyclic  stress  ampli- 
tudes, stress  ratios  (R  = minimum  stress/maximum  stress),  temperatures  and  cyclic  dwell 
times  and  stress  levels.  Nonlinear  finite  element  analyses  defined  the  disk  and  specimen  notch 
response  to  the  nominal  loading  conditions.  The  finite  element  analyses  were  used  to  cali- 
brate a simple  notch  model  actually  used  in  the  life  prediction  method.  Resulting  values  of 
notch  strain  range  and  mean  stress  were  used  to  characterize  LCF  life  exhaustion. 

A comprehensive  LCF  life  exhaustion  computer  program  was  written  utilizing  the  resulting 
notch  stress/strain  response  mooel,  the  LCF  life  exhaustion  model  including  a cumulative 
damage  effects  algorithm,  and  the  dwell  effects  model  that  were  developed  by  the  research 
program.  The  computer  program  estimates,  on  a cold  (fan)  or  hot  (turbine)  component 
basis,  LCF  life  exhaustion  as  a function  of  component  nominal  mission  load  history,  notch 
configuration,  and  material  properties. 


1.2  PROGRAM  RELEVANCE 

One  of  the  major  benefits  possible  through  the  development  of  a comprehensive  tool  for 
predicting  component  life  exhaustion  would  be  the  improvement  in  premature  retirement 
of  disks.  Presently,  the  Air  Force  is  replacing  disks  from  engines  which  show  no  evidence 
of  measurable  damage,  based  upon  existing  nondestructive  inspection  (NDI)  methods.  How- 
ever, any  method  which  promises  greater  part  life  inherently  contains  the  risk  of  premature 
failure.  Consequently,  the  results  of  this  program  need  to  be  carefully  weighted  to  ensure 
that  application  of  the  predictive  tool  does  not  increase  failure  probability  for  components 
which  currently  are  accurately  accounted  for  under  existing  systems.  Especially  important 
in  this  regard  are  high  criticality  components,  such  as  disks,  where  failure  can  cause  loss  of 
the  entire  engine. 


Application  of  the  prediction  tool  promises  a better  understanding  of  the  application  of 
testing  procedures  to  verify  accumulated  damage  in  real  hardware.  Accelerated  test  proce- 
dures which  incorporate  component  normal  operational  fleet  life  exhaustion  events  on  a one- 
to-one  basis  are  impossible  to  design.  A realistic  test  which  verifies  design  goals  requires  a 
knowledge  of  how  accelerated  endurance  damage  accumulation  is  related  to  real  time  fleet 
damage  accumulation.  Procedures  developed  within  this  program  could  result  in  significant 
cost  savings  in  designing  endurance  life  verification  tests.  Envelopment  of  these  test  proce- 
dures promises  savings  in  the  elimination  of  redesigns  during  the  development  phase  of  engine 
programs,  as  well  as  savings  in  avoidance  of  costly  redesigns. 

Recent  occurrences  of  component  failures  resulting  from  fatigue  related  damage  have  had  an 
impact  on  engine  development  programs  and  fleet  deployment.  An  example  is  the  disk  lug 
ruptures  which  occurred  during  MQT  endurance  testing  of  the  FI  00  engine.  The  MQT 
testing  was  subsequently  judged  to  be  very  severe  when  compared  to  actual  field  usage.  In 
other  instances,  cracking  of  engine  components  has  occurred  within  a few  hundred  hours  of 
operational  use,  even  though  thousands  of  hours  of  test  stand  running  were  accumulated  on 
the  engine  model.  The  cost  of  retrofit  of  component  hardware  following  service  introduction 
and  the  attendant  delays  in  full  deployment  are  major  concerns  to  both  the  services  and  the 
manufacturer. 


No  less  important  is  the  long  range  impact  of  replacement  of  life  exhausted  components.  In 
the  last  20  years,  the  cost  of  military  fighter  engines  has  multiplied  by  a factor  of  three  to 
five,  based  on  adjustment  to  1978  dollars. 

Costs  of  replacing  disks  have  risen  accordingly,  especially  for  complex  turbine  disks.  These 
increases  are  dictated  largely  by  the  requirement  for  higher  stresses  and  lighter  weight;  this 
requirement  necessitates  use  of  advanced  alloys  and  machining  of  complex  shapes.  It  can  be 
extrapolated  from  Figure  2 that  replacement  costs  of  future  engines  could  reach  five  to  ten 
times  today’s  fleet  replacement  costs.  This  program  provides  the  essential  effort  to  obtain  a 
comprehensive  life  exhaustion  diagnostic  tool  to  minimize  the  impact  of  the  escalating  costs 
of  component  replacements. 


Year 


Figure  2 Increased  Replacement  Costs  of  Compressor  and  Turbine  Disks  Are  Due  Largely 
to  Requirements  for  Higher  Stressed,  Lighter.  More  Complex  Disks 
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1 .3  MAJOR  PROGRAM  ACCOMPLISHMENTS 


The  Structural  Life  Prediction  and  Analysis  Technology  research  program  has  successfully 
established  a basis  for  LCF  life  exhaustion  modeling  relevant  to  Air  Force  needs.  It  was 
established  that  unnotched  SC  specimens,  together  with  a model  for  a workhardened  surface 
layer  in  a machined  component,  can  be  used  to  predict  the  LCF  life  of  machined  BH 
specimens.  Further,  nonlinear  cumulative  damage  effects  were  clearly  established  for  misr 
sion-relevant  load  histories,  and  these  effects  have  been  successfully  modeled  using  a fracture 
mechanics  based  methodology.  The  effects  of  dwell  (creep)  on  LCF  life  were  predicted  using 
a simple  mean  cyclic  stress  relaxation  model.  These  accomplishments  form  the  basis  of  the 
cyclic  life  prediction  algorithm. 

Cost-effective  numerical  modeling  of  the  bolt-hole  notch  was  developed  by  calibration  of 
a simple  notch  model  with  detailed,  nonlinear  finite  element  results.  The  notch  model 
algorithm  defines  the  cyclic  stress  (strain)  conditions  in  terms  of  the  strain  range  and  mean 
stress  for  the  material  adjacent  to  the  notch.  The  research  program  clearly  established  the 
validity  of  these  variables  for  predicting  LCF  life  exhaustion. 

Finally,  the  resulting  LCF  life  exhaustion  model,  incorporating  all  of  these  advanced  fea- 
tures, successfully  correlated  the  full  scale  component  test  results.  Mission  simulation  testing 
of  both  sets  of  engine  disks  included  important  variations  of  load  level,  mission  complexity, 
and  mission  ordering.  The  resulting  model  has  been  converted  to  a computer  program  for 
Air  Force  use. 

The  following  report  summarizes  these  accomplishments  and  describes  the  analytical  and 
empirical  models  in  detail.  The  computer  program  is  described  in  separate  documents*'^  *. 


*Superscript  numbers  are  references  (see  page  109). 
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SECTION  II 


MISSION  SURVEYS  AND  GENERIC  MISSIONS 

In  order  to  define  specimen  and  component  test  parameters  appropriate  to  the  investigation 
of  bolt-hole  fatigue  behavior  as  experienced  by  service  components,  an  extensive  survey  of 
missions  typically  flown  by  F 1 1 1 and  F 1 5 aircraft  was  undertaken.  The  results  of  these  sur- 
veys and  the  selection  of  the  relevant  loading  cycles  for  each  phase  of  the  effort  are  discussed 
in  this  section. 

2. 1 PHASE  I - FAN  DISK  MISSIONS 

Phase  I of  the  contract  effort  is  restricted  to  the  fan  disk  bolt-hole  region.  This  restriction 
permits  us  to  neglect  thermally  induced  stresses  and  relate  bolt-hole  loading  and  engine 
operating  conditions  solely  in  terms  of  low  rotor  speed,  N j . On  this  basis,  a typical  mission 
can  be  readily  divided  into  distinct,  activity  associated  blocks  of  subcycle  activity,  e.g.,  the 
‘Touch  and  Go”  activity  of  Figure  3. 

Table  1 lists  specific  subcycle  activities  identified  and  characterizes  each  in  terms  of  subcycle 
speed  excursion  and  frequency  of  occurrence  for  each  aircraft.  The  data  are  derived  from  a 
number  of  usage  surveys  performed  by  the  Air  Force  and  by  Pratt  & Whitney  Aircraft  and 
is  judged  to  reflect  the  mix  of  subcycle  activity  in  the  overall  Air  Force  Fleet.  Subcycle  am- 
plitudes and  levels  are  identified  as  percentage  of  take-off  maximum  low-pressure  rotor 
speed  (N  j).  The  representation  permits  consistent  comparison  of  subcycles  experienced  by 
different  disks  from  different  engines.  The  objective  was  to  set  ranges  of  stress  variation 
(proportional  to  Nj^  ) which,  when  applied  to  the  test  specimens  and  disks,  reflect  the  range 
of  subcycle  bolt-hole  loading  activity  experienced  in  the  Air  Force  Fleet. 

Those  subcycles  which  occur  at  least  100  times  during  100  typical  fleet  missions  are  shown 
in  Figure  4.  From  this  figure  it  is  seen  that  the  great  majority  (over  90  percent)  of  speed  ex- 
cursions have  a maximum  speed  very  nearly  equal  to  the  maximum  take-off  speed.  Thus,  an 
‘‘operating  line”  is  defined  along  which  nearly  all  significant  subcycle  activities  lie. 

Three  generic  mission  profiles.  Figure  5,  have  been  defined  utilizing  the  subcycle  definition 
data  discussed  above.  These  profiles  do  not  represent  usage  of  a specific  aircraft  performing 
identified  activities  but  are  instead  a composite  usage  profile  which  reflects  the  range  of  sub- 
cycle level  and  frequency  in  the  Air  Force  Fleet.  This  connectivity  is  made  evident  in  Figure  4 
where  the  generic  mission  subcycles  can  be  compared  to  the  activities  experienced  in  actual 
operation.  In  Figure  6,  the  missions  are  also  compared  to  an  operating  band  developed  by  a 
Monte  Carlo  simulation  of  Air  Force  FI  5 activity.  The  Monte  Carlo  simulation  considers 
only  major  subcycle  activity  and  corresponds  very  nearly  with  our  Mission  B.  Generic  Mission  A 
shows  a higher  subcycle  density  since  it  includes  a large  number  of  smaller  excursions  which 
the  simulation  neglects. 
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MACH  NUMBER 


TABLE  1 


SUBCYCLE  ACTIVITY  CHARACTERIZATION 

Excursion 

Excursion 

Number  of 

Max  Nj 

Min  Nj 

Excursions 

Aircraft  Engine 

(%ofT/ONi) 

(%ofT/ONj) 

per  1 00  Flights 

Fin  TF3(FP-7 

100 

0 

167 

98 

58 

100 

98 

43 

466 

85 

60 

100 

1 

1 

85 

55 

133 

T 

TF30-P-3 

100 

0 

100 

100 

74 

2825 

98 

45 

475 

97 

53 

100 

1 

I 

97 

74 

100 

f 

TF30-P-9 

100 

0 

100 

97 

62 

2825 

95 

44 

475 

98 

53 

100 

91 

65 

100 

TF- 

50-P-100 

100 

0 

100 

96 

52 

63 

88 

52 

63 

96 

88 

450 

95 

86 

450 

98 

59 

360 

75-90 

75-90  Random 

* 

97 

67 

62 

93 

82 

16 

97 

89 

10 

78 

62 

10 

100 

78 

15 

97 

62 

360 

F- 

15  Fl( 

K) 

100 

0 

100 

100 

49 

100 

99 

41 

200 

100 

50 

700 

•This  is  the  random  portion  of  TFR  Activity,  and  therefore  it  is  difficult  to  assign  a number 
of  significant  excursions  to  the  significant  N j variations. 


7 


ALTERNATING  ROTOR  SPEED 
1%  OF  TAKE-OFF  ROTOR  SPEED) 


□ SUBCYCLES  OCCURRING  AT  LEAST 
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g TFR  ACTIVITY 

generic  mission  SUBCYCLES 
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85  90  95 

MEAN  ROTOR  SPEED  1%  OF  TAKE  OFF  ROTOR  SPEED) 

Figure  4 Subcycles  Which  Include  Over  90  Percent  of  Speed  Excursion  Define  an 
“Operating  Life" 
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.TAKEOFF  CYCLE  B1 


MISSION  B 


MISSION  C 


Figure  5 Generic  Mission  Loading  Spectra  for  Test  Definition 


Figure  6 Comparison  of  Generic  Missions  and  Simulated  Usage  in  Terms  of  Accumulated 
Subcycles  Versus  Number  of  Flights 


The  loading  spectra  of  Figure  5 have  been  used  as  the  basis  for  the  test  program  including 
both  specimens  and  components.  In  addition  a computer  generated  random  loading  sequence 
corresponding  to  a ±7  percent  power  lever  angle  (PLA)  excursion  was  used  to  investigate 
damage  accumulation  during  a terrain  following  radar  (TFR)  activity.  Specific  operating 
line  levels  were  chosen  to  provide  for  loading  levels  and  test  durations  of  interest  and  will 
be  detailed  in  subsequent  sections  of  this  report. 

2.2  PHASE  II  - TURBINE  DISK  MISSIONS 

Turbine  disk  bolt-hole  stresses  respond  to  both  mechanical  loading  (high-pressure  rotor 
speed,  N2)  and  local  temperature  gradients.  Unlike  the  fan  disk  bolt-hole  response  discussed 
in  the  previous  section,  turbine  disk  stress  excursions  during  a mission  are  not  generally  pro- 
portional to  N2^.  In  order  to  properly  evaluate  the  mission  usage  of  turbine  disk  bolt  holes, 
a complete  flight  analysis,  including  the  effects  of  changing  rotor  speed  and  temperature 
gradients,  is  essential.  Such  analyses  have  been  performed  for  four  turbine  disks  and  serve 
as  the  basis  for  selection  of  relevant  test  loadings  in  this  program. 

Figure  7 shows  the  stress  response  of  the  TF30-P-100  third-stage  turbine  disk  bolt  hole  as 
used  in  the  Training  Mission  shown  in  Figure  3.  Also  noted  in  Figure  7 are  the  approximate 
bolt-hole  operating  temperatures  during  the  various  mission  activities.  The  effect  of  changing 
disk  thermal  gradients  during  an  activity  are  apparent  in  the  “bomb-run”  and  “touch  and 


Nominal  B.H.  stress,  mi  Power  setting 
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go”  activities.  Rapid  air  maneuvers  can  set  up  transient  thermal  gradients  which  may  require 
several  cycles  to  reach  a steady  state  condition.  It  is  also  apparent  that  simple  scaling  by 
N2^  would  not  predict  the  actual  stress  response.  The  principal  events  in  the  stress  history 
can  be  grouped  as  follows:  1)  a major  cycle  consisting  of  a take-off  and  landing;  2)  subcycles, 
of  varying  severity  (for  example,  the  touch  and  go  activity  constitutes  a major  subcycle 
while  the  TFR  activity  results  in  a minor  subcycle);  3)  periods  of  “dwell”  at  high  stress 
occurring  during  “climb”  activities;  and  4)  periods  of  “dweU”  at  low  stress  occurring  during 
“approach”  or  “taxi”  activities. 


Figure  7 TF30-P-1 00  Low-Pressure  Turbine  Disk  Stress  History  During  a Training  Mission 


I Principal  stress  events  for  the  four  disks  analyzed  are  shown  in  Table  2.  The  TFR  activity  is 

' j listed  as  the  minor  subcycle  for  the  TF30-P-100  low-pressure  turbine  (LPT)  disk  because  its 

* amplitude  is  comparable  to  both  TFR  and  bomb-run  activities,  thus  representing  a significant 

I damage  event  which  is  substantially  different  from  the  major  subcycle.  It  should  also  be 

i noted  that  since  actual  TF30-P-9  and  TF30-P-1 00  high-pressure  turbine  (HPT)  disks  do  not 

) j have  bolt  holes,  the  nominal  stresses  and  temperatures  listed  in  Table  2 were  computed  at 

a typical  bolt  circle  radius. 

||  Three  generic  missions  were  chosen  to  represent  the  principal  stress  events  listed  in  Table  2. 

j These  missions  are  shown  in  Figure  8,  and  their  principal  events  for  the  Ferris  Wheel  (F/W) 

disk  used  in  this  program  are  also  listed  in  Table  2.  Absolute  levels  of  stress  will  be  discussed 
in  subsequent  sections  of  this  report.  Low  stress  dwell  is  not  included  in  these  missions 
because  the  low  stress  and  temperature  levels  were  judged  insufficient  to  produce  significant 
time  dependent  material  response.  The  generic  missions  serve  to  investigate  the  interaction 
of  dwell  periods  and  subcycles  and  the  effect  of  stress  “overloads”  on  the  dwell  effect. 

II 

i|  : 
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PRINCIPAL  STRESS  EVENTS  OF  TURBINE  DISKS  AND  GENERIC  MISSIONS 


Mission 


SECTION  III 


TEST  PROGRAM  OVERVIEW 

The  development  of  the  cumulative  damage  models  in  the  contract  involved  three  major 
areas  of  fatigue  testing  and  prediction.  First,  basic  material  fatigue  data  were  obtained  for 
simple  cycle  loading.  Cumulative  damage  specimen  tests  were  conducted  next  to  identify 
the  cumulative  damage  mechanisms.  Finally,  the  cumulative  damage  models  were  substan- 
tiated with  complex  mission  testing  of  component-like  disks. 

3. 1 SIMPLE  CYCLE  FATIGUE  TESTS 

The  Strain  Controlled  (SC)  specimen  shown  in  Figure  9 was  used  to  develop  a basic  simple 
cycle  fatigue  life  prediction  equation.  This  specimen  was  used  because  gage  section  stresses 
and  strains  are  known  explicitly  for  each  test  and  thus  allow  the  development  of  a fatigue 
system  which  is  independent  of  notch  geometry  and  notch  stress/strain  estimation.  The  SC 
specimens  were  electrochemically  machined  to  achieve  a stress  free  surface. 


Figure  9 Strain  Controlled  (SC)  Specimen  ( 77-551-9091 ) 
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The  Bolt  Hole  (BH)  specimens  shown  in  Figures  10  and  1 1 were  used  to  find  the  relation 
between  the  simple  cycle  LCF  behavior  of  bolt  holes  and  the  smooth  SC  specimens.  From 
these  back-to-back  tests,  the  presence  of  a bolt-hole  surface  residual  stresses  was  inferred  in 
both  the  PWA  1216  and  PWA  1057  bolt  holes.  A surface  residual  stress  algorithm  which  ac- 
counts for  the  effect  of  the  residual  stress  and  the  cyclic  hardened  condition  of  the  surface 
layer  was  developed  and  is  described  in  Appendix  B.  The  use  of  this  algorithm  showed  that 
the  basic  fatigue  behavior  equation  developed  with  the  SC  specimen  was  applicable  to  the 
BH  specimens  as  well.  This  merging  of  the  SC  and  BH  fatigue  data  confirmed  the  expecta- 
tion that  fatigue  behavior  prediction  can  effectively  be  transferred  from  specimens  to  com- 
ponents if  the  local  notch  mechanical  conditions  are  accurately  known. 


5.50  NOM, 


1.446  NOM, 


Figure  10  Phase  / Bolt  Hole  ( BID  Spec  imen,  Kj 


Figure  11  Phase  II  Bolt  Hole(BH)  Specimen,  Kj  = 2.34  ( XPN-150211 


The  simple  cycle  bolt-hole  test  conditions  were  chosen  to  simulate  principal  component 
usage  identified  in  the  mission  surveys  discussed  in  Section  II.  Bolt-hole  specimen  loading 
cycles  are  shown  schematically  in  Figure  12. 

In  Phase  I testing,  the  principal  features  were  major  cycle  tests  corresponding  to  take-off 
and  landing  and  various  amplitudpc  of  subcycle  testing  corresponding  to  in-flight  activities. 
Two  levels  of  maximum  nominal  stress  evaluated  in  the  major  cycle  tests,  and 

two  levels  of  nominal  stress  rar.ge,  AS,  were  tested  in  the  subcycle  tests.  All  tests  were  con- 
ducted at  room  temperature. 

In  Phase  II,  major  cycle  tests  with  and  without  the  effect  of  high  stress  dwell,  and  subcycle 
tests  without  dwell  were  conducted.  Three  levels  of  S^^j^  and  three  levels  of  AS  were  eval- 
uated. All  tests  were  conducted  isothermally  at  900° F to  simulate  engine  environment. 

Details  of  the  testing  and  fatigue  model  are  discussed  in  Section  V. 
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Major  cycle  tests 


Subcycle  tests 


Phase  I 


During  Phase  I,  a limited  number  of  simple  cycle  LCF  tests  were  conducted  using  the  notched 
round  bar  (NRB)  specimen  shown  in  Figure  13.  However,  elastoplastic  finite  element  analysis 
of  this  specimen  showed  that  the  stress-strain  response  of  the  notch  was  significantly  different 
in  character  from  the  response  of  either  the  SC  or  BH  specimens  due  to  a very  high  degree  of 
biaxiality  in  the  notch.  These  results  were  previously  reported^.  The  incorporation  of  the 
data  W3s  therefore  inappropriate  for  the  development  of  a bolt-hole  life  prediction  model. 
Consequently,  no  NRB  tests  were  conducted  during  Phase  II  of  this  contract,  and  NRB 
behavior  will  not  be  considered  in  this  report. 


3.2  CUMULATIVE  DAMAGE  TESTS 

Cumulative  damage  investigations  were  conducted  using  various  combinations  of  test  condi- 
tions already  characterized  in  the  simple  cycle  testing  program.  Only  two  sorts  of  cumulative 
damage  tests  are  possible;  block  loading  or  sequenced  tests  as  shown  schematically  in  Figure 
14.  Nearly  all  testing  was  of  the  block  loading  type  in  which  load  c rdering  effects  can  easily  be 
observed.  The  sequenced  test  does  not  evaluate  the  effect  of  order  ng. 
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Figure  13  Notched  Round  Bar(NRB)  Specimen 


Figure  14  Cumulative  Damage  Bolt  Hole  Tests 


In  both  Phase  I and  Phase  II,  the  block  loading  tests  revealed  a strong  influence  due  to  order- 
ing which  was  modeled  successfully  by  double-damage  concepts.  In  Phase  II,  the  effect  of 
the  dwell  cycle  was  also  found  to  be  modeled  properly  by  double-damage  concepts. 

The  test  conditions  were  chosen  to  reflect  the  ordering  effects  of  interest  identified  by  the 
mission  surveys.  For  example,  the  effect  of  a stress  overload  following  a dwell  cycle  was  in- 
vestigated in  both  block  loading  and  sequenced  tests  during  Phase  II. 

Details  of  the  cumulative  damage  specimen  testing  are  discussed  in  Section  VI. 
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3.3  FERRIS  WHEEL  DISK  TESTS 


The  cumulative  damage  models  developed  in  the  specimen  test  program  were  verified  by  test- 
ing of  full-size  disk-like  components  in  a “Ferris  Wheel”  (F/W)  which  simulates  engine  loading 
through  radially  applied  loads  at  the  disk  rim.  The  Phase  I tests  were  conducted  at  room  tem- 
perature, while  the  Phase  II  tests  were  conducted  at  900°  F without  temperature  gradients 
at  the  bolt  circle. 

Four  disks  were  tested  in  Phase  I and  three  disks  in  Phase  II.  The  generic  missions  shown  in 
Figures  5 and  8 were  used  as  the  basis  for  disk  loading.  The  effects  of  mission  mixture  and 
different  maximum  load  levels  were  included  to  provide  a reasonably  comprehensive  test  of 
the  cumulative  damage  prediction  model. 

The  F/W  tests  are  discussed  in  detail  in  Section  VI. 
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SECTION  IV 


CHARACTERIZATION  OF  LOCAL  NOTCH  BEHAVIOR 

Accurate  characterization  of  local  bolt-hole  mechanical  conditions  such  as  strain  range  and 
mean  stress  is  fundamental  to  accurate  LCF  life  prediction  of  the  hole.  Because  the  oper- 
ating load  levels  for  engine  disks  often  result  in  local  plastic  deformation  at  the  bolt  hole, 
elastoplastic  stress  analyses  are  essential.  In  addition,  the  operating  temperatures  and  stresses 
in  turbine  disk  bolt  holes  are  often  sufficiently  high  to  cause  time  dependent  material  non- 
linear behavior.  Calculation  of  local  notch  stresses  and  strains  in  the  presence  of  these  two 
nonlinear  effects  are  described  below. 

4.1  ELASTOPLASTIC  ANALYSIS 

Two  nonlinear  models  were  used  in  this  contract  to  predict  local  mechanical  variables  for 
both  monotonic  and  cyclic  plastic  response.  The  first  model  is  the  well-known  Neuber 
model^ . The  stress  concentration  and  strain  concentration  factors,  given  by 

K^  = a/S, 

Kg  = eE/S 

are  related  by  the  elastic  value  of  the  stress  concentration  factor, 

K„Kg  = Kt2  (1) 

Equation  (1 ) defines  a Neuber  hyperbola  for  a given  applied  nominal  stress  (S).  The  inter- 
section of  the  Neuber  hyperbola  and  the  stress-strain  curve  defines  the  local  notch  stress  and 
strain  values  as  shown  in  Figure  1 5. 


Figure  15  Determination  of  Local  Cyclic  Conditions  by  the  Neuber  Calculation 
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The  second  model  is  an  elastoplastic  finite  element  analysis  using  the  MARC*  computer 
program.  The  finite  element  analysis  was  judged  to  be  the  most  accurate  characterization 
method,  but  was  too  cumbersome  for  routine  calculations.  The  Neuber  notch  model  was 
judged  to  be  extremely  efficient,  but  had  to  be  modified  by  the  finite  element  results  to 
obtain  necessary  levels  of  modeling  accuracy. 

Comparison  of  the  two  models  for  the  BH  specimens  is  shown  in  Figures  16  and  17  in  terms 
of  the  local  strain  concentration  (Kg)  and  stress  concentration  (K^j)  factors.  The  comparison 
is  made  at  a material  location  just  below  the  bolt-hole  surface  which  corresponds  to  file  loca- 
tion of  the  nearest  finite  element  stress  solution  point.  The  comparison  is  judged  to  be  valid 
also  for  the  bolt-hole  surface  itself.  The  comparison  is  given  for  both  the  monotonic  stress- 
strain  curve  and  the  hysteresis  curve  for  the  material  tested  in  Phase  I and  for  the  monotonic 
curve  used  in  Phase  II.  Details  of  the  stress-strain  curves  are  given  in  Appendix  A.  It  is  prop- 
erly concluded  that  the  Neuber  model  gives  good  results  for  local  stress  and  strain  values  so 
long  as  the  stress-strain  curve  exhibits  a high  degree  of  strain  hardening  (i.e.,  the  stress-strain 
curve  is  relatively  “steep”)  and  can,  therefore  be  used  with  the  hysteresis  curves.  However, 
use  of  the  Neuber  model  with  the  monotonic  curves  (low  strain  hardening)  is  seen  to  lead  to 
significant  over-estimation  of  strains.  Figures  16  and  17  can  be  used  to  modify  or  correct 
the  Neuber  calculation.  The  Neuber  correction  curves  for  the  bolt-hole  specimens  responding 
to  the  monotonic  stress-strain  curves  are  given  in  Figure  18.  The  Neuber  approach,  properly 
corrected  for  use  with  the  monotonic  stress-strain  curves,  was  used  in  both  Phase  I and 
Phase  II  of  the  contract. 


Figure  16  Finite  Element  Results  Versus  Neuber  Calculation  for  PWA  1216  Materidl 


*The  MARC  general  purpose,  nonlinear  finite  element  code,  MARC  Analysis  Corporation, 
Providence,  R.  I. 
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Figure  1 7 Finite  Element  Results  Versus  Neuber  Calculation  for  PWA  1057  Material 


Strain  from  Neuber  calculation  (n/n.) 


Figure  18  Strain  Calibration  Curves  for  Neuber  Calculation 
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4.2  BOLT-HOLE  SURFACE  LAYER  RESPONSE 

In  both  Phase  I and  Phase  II  of  the  contract,  recognition  of  a bolt-hole  surface  layer  was 
crucial.  Machining  of  the  BH  specimens  and  F/W  disk  bolt  holes  is  known  to  introduce  com- 
pressive surface  residual  stresses.  A surface  layer  response  algorithm  was  developed  which 
accounts  for  the  effects  of  an  initial  stress  state,  a^,  and  the  cyclic  hardened  condition  de- 
veloped in  the  surface  layer  during  the  machining  operation.  Appendix  B provides  an  explan 
ation  of  this  algorithm. 

The  modified  Neuber  approach  described  in  the  previous  section  was  used  in  conjunction 
with  this  surface  layer  algorithm  to  calculate  the  stress  response  of  the  surface  layer  and  the 
underlying  subsurface  material  during  fatigue  testing.  Figure  19  illustrates  the  method  for 
two  cycles  of  loading  of  a bolt-hole  specimen.  Stress  and  strain  values  during  the  cycle  for 
this  example  are  listed  in  Table  3.  The  modified  Neuber  approach  is  used  to  calculate  the 
subsurface  stress  and  strain  at  point  A since  the  low-hardened,  monotonic  curve  is  involved. 
SubscQuent  stress  excursions,  from  A to  B and  from  B to  C for  both  surface  and  subsurface, 
follow  the  high-hardened,  hysteresis  curve,  so  the  Neuber  calculation  need  not  be  modified. 


Load  history  Subsurface  response 


Figure  1 9 Surface  Layer  Plastic  Response 


TABLE  3 


SURFACE  AND  SUBSURFACE  STRESSES  AND  STRAINS 
FOR  ILLUSTRATIVE  EXAMPLE 


Surface  Layer 

Reference  Point  — " 

in  Load  History  Stress  Strain 


Subsurface  Layer 


Stress  Strain 


O 

A 

B 

C 


‘^o  ‘o  O 

+ e + e^  S(i) 

Og+a  (€)-o  (Ae)  ?-Ae  + e a(^-o(Ac) 

a^+a(^)  i + e d 


O 


Ae 


4.3  STRESS  RELAXATION 


The  temperatures  and  stress  levels  experienced  in  turbine  disks  and  in  the  test  conditions  of 
Phase  II  were  sufficiently  severe  to  cause  noticeable  creep  in  a standard  creep  specimen  as 
reported  in  Appendix  A.  In  a strain  controlled  environment,  the  material’s  nonlinear  behavior 
will  be  exhibited  as  a stress  relaxation.  Two  cyclic  relaxation  tests  of  the  Phase  II  material, 
PWA  1057,  were  conducted  using  the  SC  specimen  shown  in  Figure  9,  These  data,  reported 
in  Appendix  A,  showed  that  the  mean  stress  during  cycling  decayed  exponentially  and  the 
relaxation  of  the  mean  stress,  could  be  written  as 


®R  = C 


(C-0002T.,) 


(2) 


where:  is  the  time  averaged  mean  stress  in  the  first  cycle  and  T is  the  total  time  spent 

at  the  maximum  strain  limit  during  cycling. 


Since  all  testing  under  Phase  II  was  isothermal,  the  assumption  is  made  that  the  exponential 
decay  rate  is  the  same  for  all  stress  levels  and  the  relaxation  is  given  by  (2).  A complete 
derivation  of  for  complex  missions  and  simple  cycle  tests  is  provided  in  Appendix  C. 

The  relaxation  of  both  surface  and  subsurface  stresses  is  recognized  fully  in  the  simple  cycle 
and  complex  cycle  LCF  life  predictions  in  the  Phase  II  effort. 
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SECTION  V 


SIMPLE  CYCLE  LIFE  MODEL 


A mathematically  simple  model  has  been  used  to  correlate  the  simple  cycle  LCF  data.  The 
model  recognizes  that  LCF  life  is  dependent  upon  the  severity  of  the  hysteresis  activity  ex- 
perienced locally,  measured  in  terms  of  the  local  strain  range,  Ae,  and  the  mean  stress  value, 
®m’  which  this  activity  occurs.  The  life  correlating  model  is  given  by 


N = AAf® 


(3) 


where  N is  the  LCF  life,  and  A,  B,  and  C are  material  dependent  constants.  Constant  life 
curves  from  (3)  are  shown  in  Figure  20. 

This  section  discusses  the  specific  use  of  (3)  in  correlating  the  simple  cycle  LCF  data  in 
Phase  I and  Phase  II. 


,|  Figure  20  Constant  Life  Curves 
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S.  1 PHASE  I - SIMPLE  CYCLE  LIFE  CORRELATION 

During  Phase  I PWA  1216  strain  controlled  specimens  (Figure  9)  of  a single  heat  code  (CAAZ) 
were  cycled  between  fixed  strain  limits  and  periodically  inspected  to  establish  the  number  of 
cycles,  N,  needed  to  initiate  a 1/32  inch  surface  length  crack.  Table  4 provides  a list  of  test 
conditions  and  resulting  LCF  lives  for  each  specimen.  Also  included  in  Table  4 are  the  cor- 
responding values  of  stress  and  strain  occurring  during  the  first  loading  and  during  subsequent 
cycling  between  the  strain  limits.  The  values  of  stress  listed  here  were  obtained  directly  from 
the  monotonic  and  hysteresis  curves  given  in  Appendix  A and  do  not  involve  any  nonlinear 
analysis.  The  mean  stress  level,  o,„,  is  given  by 

On,  = o - (Ae)  (4) 


where  o is  the  maximum  stress  occurring  during  the  first  loading  cycle  obtained  from  the 
monotonic  stress  strain  curve,  and  a(Ae)  is  the  cyclic  stress  resulting  from  the  strain  excur- 
sion Ae  obtained  from  the  hysteresis  stress  strain  curve. 
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The  logarithmic  form  of  (3), 

log  N = log  A + B log  Ae  + COj^  (5) 

lends  itself  to  standard  multivariate  linear  regression  techniques  in  which  log  Ae  and  Opi 
treated  as  the  independent  variables.  Appendix  D contains  a discussion  of  multivariate  re- 
gression analyses.  The  constants  log  A,  B,  and  C determined  by  such  an  analysis  can  be  iden- 
tified in  the  following  equation. 

log  N = -8.0979  - 6.7319  log  Ae  - 0.020306  (6) 

TABLE  4 


PHASE  I STRAIN  CONTROLLED  SPECIMEN 
SIMPLE  CYCLE  LCF  DATA 


Strain  Limits 

First ! 

Loadine 

Subsequent  Loading 

Actual 

Predicted 

Test 

^min 

^max 

t 

Ae 

a(Ae) 

‘'m 

Life 

Life* 

Condition 

(in/in) 

(in/in) 

(ksi) 

(inyin) 

(in/in) 

(ksq 

(ksO 

(cycles) 

(cycles) 

SCI 

0.0000 

0.0157 

179.9 

0.0157 

0.0157 

253.3 

53.3 

740 

960 

0.0000 

0.0157 

179.9 

0.0157 

0.0157 

253.3 

53.3 

750 

960 

0.0000 

0.0157 

179.9 

0.0157 

0.0157 

253.3 

53.3 

850 

960 

SC2 

0.0030 

0.0155 

179.8 

0.0155 

0.0125 

212.0 

73.8 

1,800 

1,700 

0.0030 

0.0155 

179.8 

0.0155 

0.0125 

212.0 

73.8 

2,000 

1,700 

0.0030 

0.0155 

179.8 

0.0155 

0.0125 

212.0 

73.8 

2,400 

1,700 

SC3 

0.0062 

0.0155 

179.8 

0.0155 

0.0093 

162.8 

98.4 

2,800 

3,800 

0.0062 

0.0155 

179.8 

0.0155 

0.0093 

162.8 

98.4 

2,800 

3,800 

0.0062 

0.0155 

179.8 

0.0155 

0.0093 

162.8 

98.4 

3,400 

3,800 

0.0062 

0.0155 

179.8 

0.0155 

0.0093 

162.8 

98.4 

3,600 

3,800 

SC4 

0.0010 

0.0103 

173.3 

0.0103 

0.0093 

162.8 

91.9 

5,700 

5,200 

0.0010 

0.0103 

173.3 

0.0103 

0.0093 

162.8 

91.9 

8,700 

5,200 

0.0010 

0.0103 

173.3 

0.0103 

0.0093 

162.8 

91.9 

9,400 

5,200 

SC5 

0.0000 

0.0093 

162.8 

0.0093 

0.0093 

162.8 

81.4 

6,600 

8,600 

0.0000 

0.0093 

162.8 

0.0093 

0.0093 

162.8 

81.4 

8,000 

8,600 

0.0000 

0.0093 

162.8 

0.0093 

0.0093 

162.8 

81.4 

11,000 

8,600 

0.0000 

0.0093 

162.8 

0.0093 

0.0093 

162.8 

81.4 

11,500 

8,600 

SC6 

0.0002 

0.0093 

162.8 

0.0093 

0.0073 

127.8 

98.9 

12,500 

19,100 

0.0002 

0.0093 

162.8 

0.0093 

0.0073 

127.8 

98.9 

18,100 

19,100 

* using  equation  (S)  with  fatigue  constants  (10) 
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The  SC  data  cortelation  and  the  95  percent  confidence  band  is  shown  in  Figure  21.  The  con- 
fidence band  shows  the  limits  of  lives  within  which  a data  point  would  be  expected  to  fall 
95  percent  of  the  time  based  on  the  accuracy  of  the  prediction  model  (6).  Figure  21  shows 
that  the  good  correlation  has  been  achieved  for  the  data  set  taken  as  a whole  and  that  the 
correlation  of  individual  groups  of  data  is  quite  consistent. 


Figure  21  Phase  I Strain  Controlled  ( SC)  Data  Correlation 

Figure  22  compares  the  Weibull  distributions  for  the  data  as  correlated  with  (6)  and  for  a 
“perfect”  correlation  model.  (See  Appendix  D for  further  discussions  of  the  Weibull  distri- 
bution .)  The  additional  scatter  introduced  by  the  correlation  model  (indicated  by  the  de- 
crease in  slope  0)  is  not  judged  to  be  excessive,  especially  in  view  of  the  relative  simplicity 
of  the  correlation  model.  Furthermore,  the  Weibull  distribution  is  continuous,  and  Table  4 
reveals  no  patterns  of  consistently  poor  correlation  for  any  particular  subset  of  the  data. 

The  overall  correlation  is  judged  to  be  very  good. 

Bolt-hole  specimens,  as  shown  in  Figure  10,  were  loaded  cyclicly  between  fixed  nominal 
stress  levels  and  periodically  inspected  to  establish  the  1/32  inch  crack  life.  Table  5 lists 
the  test  conditions  and  actual  lives.  Included  in  Table  5 are  the  stresses  and  strains  used 
to  characterize  the  local  notch  behavior.  The  values  are  calculated  using  the  appropriate 
modified  Neuber  calculation  and  the  surface  residual  stress  algorithms  discussed  in  Section 
IV  and  Appendix  B.  Figure  23  illustrates  the  hysteresis  activities  for  each  test  condition  in 
a-€  space.  The  mean  stress  of  the  surface  layer  is  given  by 

= Oo  + a (?)  - -^  a (Ae)  (7) 

The  value  of  is  not  known  explicitly  but  was  found  through  regression  techniques  as  fol- 
lows. 


The  mean  stress  for  both  SC  and  BH  specimens  is  composed  of  two  parts;  an  initial 
stress  state,  and  o^,  the  stress  excursions  from  the  initial  stress  state  to  the  mean  stress  level. 
That  is 


'm 


= + o. 


(8) 


where 


1 

2 


a (Ae) 


a 


a 


o(e)  - 


1 

2 


a (Ae) 


for  SC  specimens 


for  BH  specimens 


For  the  electrochemically  machined  SC  specimens,  Oq  - 0 while  0 but  is  unknown  for 
the  BH  specimens.  (5)  can  be  rewritten  as 

log  N = log  A + B log  Ae  + C (D  Oq  + O3)  (9) 

where  D = 0 for  SC  specimens 
D = 1 for  BH  specimens 

This  equation  reduces  to  (3)  for  the  SC  specimens.  A multivariate  regression  analysis  with 
three  independent  variables,  (log  Ae,  a^,  and  D)  for  each  data  point,  will  yield  the  constants 
log  A,  B,  C and  Ca^  for  the  combined  BH  plus  SC  data  set.  The  value  of  can  then  be 
easily  calculated.  The  values  for  such  a regression  analysis  are  as  follows; 

log  A = -8.0508 

B = -6.7195  (10) 

C = -0.020476 

Oq  = -136.1 

It  is  clear  that  the  material  fatigue  constants,  log  A,  B,  and  C,  agree  very  closely  with  those 
obtained  from  considering  the  SC  data  only,  (6).  Predicted  values  of  life  using  constants 
from  the  combined  data  set,  are  shown  in  Tables  4 and  5.  Figure  24  shows  that  the  overall 
correlation  of  the  combined  data  set  is  excellent;  all  data  points  lie  well  within  the  95  per- 
cent confidence  band  expected  for  SC  data  only,  and  the  mean  lives  of  each  data  set  are  cor- 
related very  well.  The  Weibull  plot  of  Figure  25  shows  that  the  scatter  in  the  overall  correla- 
tion is  not  significantly  larger  than  the  inherent  scatter  in  the  data.  It  is,  in  fact,  closer  to 
the  scatter  in  a perfect  correlation  than  was  achieved  by  the  SC  data  alone,  see  Figure  22. 
This  is  due  to  a better  definition  of  rhe  scatter  resulting  from  a larger  population  and  the 
obvious  success  of  the  merger  of  the  BH  and  SC  data  sets.  Thus,  the  correlation  is  judged  to 
be  very  good,  and  it  is  clear  that  the  correlation  model  adequately  accounts  for  the  major 
fatigue  influencing  parameters. 


Scatter  in  Phase  I 
SC  correlation 


Scatter  in  a 
perfect  correlation 
yS  =7.0 


Actual  life 

Avg.  life  at  each  condition 


Figure  22  Phase  I SC  Correlation  Weibull  Plots 


TABLE  5 

PHASE  I BOLT-HOLE  SPECIMEN 
SIMPLE  CYCLE  LCF  DATA 


Nominal 

Stress 

First  Loading 
Subsurface  Surface 

Subsequent  Loading 

Actual 

Predicted 

Test 

Limits 

S 

o(f) 

Ae 

0 (Ae) 

"m’^o 

Life 

Life* 

Condition 

(Itsi) 

(Itsi) 

(in/in) 

(ksi) 

(in/in) 

(Itsi) 

(ksi) 

(cycles) 

(cycles) 

BHI 

50  ± 50 

181.5 

0.0191 

279.0 

0.0151 

246.1 

156.0 

4,500 

6,000 

50  ± 50 

I8I.5 

0.0191 

279.0 

0.0151 

246.1 

156  0 

5,300 

6,000 

50±  50 

I8i.5 

0.0191 

279.0 

0.0i51 

246.1 

156.0 

5,600 

6,000 

50±  50 

I8I.5 

0.0191 

279.0 

0.0151 

246.1 

156.0 

5,800 

6,000 

50±  50 

181.5 

0.0191 

279.0 

0.0151 

246.1 

156  0 

6,000 

6,000 

50  ±50 

181.5 

0.0191 

279.0 

0.0151 

246.1 

156.0 

6,500 

6,000 

50±  50 

181. 5 

0.0191 

279.0 

0.0151 

246.1 

156.0 

7,200 

6,000 

50±  50 

181.5 

0.0191 

279.0 

0.0151 

246.1 

156.0 

7,400 

6,000 

50±  50 

I8I.5 

0.0191 

279.0 

0.0151 

246.1 

156.0 

7,600 

6,000 

BK2 

60  ±40 

181.5 

0.0191 

279.0 

0.0119 

203.0 

177.5 

10,000 

10,800 

60  ±40 

181.5 

0.0191 

279.0 

0.0119 

203.0 

177.5 

13,900 

10,800 

BH3 

70  ±30 

181.5 

0.0191 

279.0 

0.0088 

157.6 

200.2 

23,500 

28,100 

70  ±30 

181.5 

0.0191 

279.0 

0.0088 

157.6 

200.2 

31,000 

28,100 

BH4 

40  ±40 

179.1 

0.0127 

214.8 

0.0119 

203.0 

113.3 

157,000 

223,000 

* using  equation  (S)  with  fatigue  constants  (10) 
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Scattgr  in  correlation 
Phase  1 combined  data  set 
/3  =5.2 


Scatter  in  a 
perfect  correlation 
/8=7.5 


Figure  25  Phase  I Combined  Data  Set  Weibull  Plots 


The  value  of  resulting  from  the  regression  must  be  evaluated  in  terms  of  its  physical  in- 
terpretation in  the  surface  layer  algorithm;  the  value  of  —136.1  ksi  is  certainly  a reasonable 
value  for  a residual  surface  stress. 

The  close  agreement  of  the  fatigue  constants  and  the  reasonable  value  of  indicate  that 
the  surface  layer  algorithm  is  reasonable.  In  order  to  further  evaluate  the  algorithm,  an  at- 
tempt was  made  to  predict  the  BH  data  using  mean  stress  values  given  by  (4)  (i.e.,  the  value 
of  Ojj^  as  if  no  surface  layer  were  present)  and  (6)  developed  for  the  SC  data  alone.  The  re- 
sults of  this  attempt  are  shown  in  Figure  26.  It  is  obvious  that  a very  significant  improve- 
ment is  achieved  by  the  use  of  the  surface  layer  residual  stress  algorithm. 

Several  important  conclusions  were  made  from  the  above  analyses: 

1)  The  local  mechanical  parameters,  strain  range,  Ae,  and  mean  stress,  as  used  in 
the  life  model  (3)  are  the  principal  variables  affecting  the  fatigue  life  in  bolt-hole- 
like regions. 

2)  The  constants  A,  B and  C are  basic  materia!  fatigue  constants  that  can  be  deter- 
mined from  well  controlled  SC  testing.  This  suggests  an  alternate  approach  to 
finding  the  value  of  Oq  for  a notch;  with  material  constants  A,  B,  and  C known 
from  SC  tests,  the  value  of  Oq  can  be  evaluated  by  using' (9).  This  approach  was 
used  in  Phase  II. 
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3)  A work  hardened  surface  layer  with  an  initial  stress  state  is  present  in  the  bolt 
holes,  and  the  surface  layer  algorithm  of  Appendix  B adequately  models  its  be- 
havior. 

4)  Accurate  calculation  of  Ae  and  is  essential  to  accurate  fatigue  prediction  and 
is  a fundamental  link  in  translating  specimen  fatigue  results  to  component  life 
prediction. 


Figure  26  Prediction  of  Phase  I BH  Data  Without  Consideration  of  Surface  Layer 

5.2  PHASE  II  - SIMPLE  CYCLE  LIFE  CORRELATION 

During  Phase  II  testing,  the  element  of  elevated  temperature  effects  was  added  to  the  inves- 
tigation of  fatigue  behavior.  All  tests  were  conducted  isothermally  at  900° F using  a single 
heat  code  (XNNZ)  of  PWA  1057.  The  basic  approach  in  the  investigation  remained  essen- 
tially the  same  as  that  followed  in  the  Phase  I room  temperature  testing.  That  is,  the  same 
life  correlating  model  (3)  was  employed,  and  SC  tests  were  used  to  establish  the  material’s 
fatigue  constants  A,  B,  and  C.  Next,  bolt-hole  testing,  with  appropriate  accounting  of  sur- 
face layer  and  stress  relaxation,  was  evaluated  against  the  fatigue  model. 
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Table  6 lists  the  individual  SC  test  conditions  and  resulting  1 /32  inch  crack  lives.  As  noted 
in  Appendix  A,  the  PWA  1057  material  hysteresis  curves  were  dependent  upon  strain  range 
and,  in  general,  showed  more  variability  than  the  PWA  1216  room  temperature  behavior, 
thus  making  calculation  of  mean  stress  levels  from  standard  material  curves  more  difficult 
than  in  Phase  I.  In  addition,  even  small  temperature  fluctuations  while  testing  in  a strain 
controlled  mode  caused  noticeable  variations  in  recorded  specimen  load  levels.  Consequently, 
the  average  mean  stress  values  for  each  test  condition  were  obtained  from  the  load  records 
of  each  test  and  used  in  the  data  reduction.  The  regression  analysis  of  these  data,  using  the 
logarithmic  form  of  (3),  yielded  the  material  fatigue  constants  as  follows: 


log  A = -6.0796 

B = -4.8057  (11) 

C » -0.013647 


TABLE  6 

PHASE  II  STRAIN  CONTROLLED  SPECIMEN 
SIMPLE  CYCLE  LCF  DATA 


Test 

Specimen 

^niin 

®max 

Ac 

Actual 

Life 

Predicted 

Life 

Condition 

Number 

(in/in) 

(in/in) 

(in/in) 

Ocsji 

(cycles) 

(cycles) 

SCI 

1057-9 

0.0000 

0.0140 

0.0140 

2.0 

900 

600 

SC2 

1057*12 

0.0050 

0.0140 

0.0090 

27.0 

2300 

2400 

1057-18 

0.0050 

0.0140 

0.0090 

41.5 

1000 

1500 

SC3 

1057-20 

0.0042 

0.0118 

0.0076 

34.5 

3000 

4300 

1057-21 

0.0041 

0.0119 

0.0078 

39.0 

1600 

3300 

1057-23 

0.0042 

0.0118 

0.0076 

42.5 

2600 

3400 

1057-25 

0.0043 

0.0121 

0.0077 

56.0 

2000 

2100 

SC4 

1057-22 

0.0000 

0.0098 

0.0098 

26.5 

2000 

1600 

1057-24 

0.0005 

0.0098 

0.0093 

42.0 

2000 

2600 

SC5 

1057-27 

0.0022 

0.0098 

0.0086 

44.0 

1000 

1800 

1057-28 

0.0022 

0.0099 

0.0087 

48.5 

2300 

1500 

SC6 

1057-36 

0.0025 

0.0090 

0.0065 

40.0 

7000 

7700 

SC7 

1057-30 

0.0003 

0.0079 

0.0076 

40.5 

6000 

3600 

1057-31 

0.0002 

0.0079 

0.0077 

32.5 

6000 

4300 

SC8 

1057-32 

0.0015 

0.0078 

0.0063 

57.0 

9000 

5200 

1057-37 

0.0015 

0.0080 

0.0065 

46.0 

8500 

6400 
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The  overall  correlation  is  shown  in  Figure  27.  The  Weibull  curves  in  Figure  28  compare  the 
degree  of  scatter  in  the  correlation  to  the  best  estimate  of  perfect  correlation.  It  is  seen  that 
the  increase  in  scatter  introduced  by  the  correlation  model  is  comparable  to  that  observed 
in  Phase  I.  And  again,  no  patterns  of  overprediction  or  underprediction  are  apparent  from 
Table  6.  As  expected,  the  higher  degree  of  variability  associated  with  high  temperature 
testing  and  material  behavior  is  reflected  in  a larger  degree  of  scatter  in  both  the  material’s 
inherent  scatter  and  in  overall  data  correlation. 

As  concluded  in  the  Phase  I testing,  the  constants  (11)  were  viewed  as  material  constants 
and  were  used  directly  in  the  correlation  of  the  BH  data.  However,  there  is  a significant  dif- 
ference in  the  nature  of  the  BH  tests  between  Phase  I and  Phase  II  which  must  be  recog- 
nized. 

The  high  temperature  environment  caused  relaxation  of  the  mean  stress  in  the  bolt-hole  re- 
gion. Since  the  mean  stress  has  been  identified  as  a primary  fatigue  correlating  variable,  re- 
laxation models  had  to  be  incorporated  in  the  characterization  of  the  mean  stress.  Of  the 
several  models  investigated,  the  model  discussed  in  Appendix  C and  in  Section  IV  was 
judged  to  be  the  most  straight-forward  and  reasonable.  This  model  was  developed  from  di- 
rect observation  of  cyclic  relaxation  tests,  and  employs  a single  decay  rate,  a (=  0.002  hr  ), 
which  is  viewed  as  characteristic  of  the  material’s  behavior  at  900°F. 
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Figure  27  Phase  II  SC  Data  Correlation 
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Figure  28  Phase  II  SC  Correlation  Weibull  Plots 


The  BH  test  conditions  were  selected  to  span  the  range  of  load  levels  and  nominal  stress 
ratios,  Rj,  pertinent  to  actual  component  usage  as  described  in  the  mission  surveys  (Section 
II).  The  test  conditions  and  actual  lives  are  given  in  Table  7,  and  the  surface  and  subsurface 
stress  conditions  are  shown  in  Figure  29. 

Prior  to  incorporating  either  the  mean  stress  relaxation  model  or  the  surface  layer  algorithm 
into  the  calculation  of  BH  mean  stresses,  an  attempt  was  made  to  predict  the  BH  data  using 
the  fatigueconstants  (1 1)  and  mean  stress  values  given  by 

= a - a (Ae)  (12) 

that  is,  as  if  neither  a surface  layer  nor  stress  relaxation  were  present  in  the  BH  specimens. 
Figure  30  shows  the  results  of  this  prediction.  It  is  clear  that  ( 1 2)  does  not  adequately  de- 
scribe the  mean  stress. 

Incorporating  the  surface  layer  mean  stress  algorithm  in  the  life  equation  for  the  BH  speci- 
mens 
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TABLE  7 


PHASE  II  BOLT  HOLE  SPECIMEN 
SIMPLE  CYCLE  LCF  DATA 


Nominal 

First  Loading 

Predicted 

Stress 

Subsurface 

Surface 

Subsequent  Loading 

Actual 

Life 

Test 

Limits 



o 



€ 

o(t) 

Ae 

a(Ae) 

ti/t2 

Life 

Eq. (16) 

Condition 

-(ksi) 

(ksi) 

(in/in) 

(ksi) 

(in/in) 

(ksi) 

(seconds) 

(cycles) 

(cycles) 

BHl 

42  ±42 

144 

0.0081 

224 

0.0068 

196 

2/2 

15,700 

55,800 

144 

0.0081 

224 

0.0068 

196 

2/2 

49,000 

55,800 

BH2 

50±  50 

148 

0.0112 

273 

0.0082 

229 

2/2 

7,000 

9,300 

148 

0.0112 

273 

0.0082 

229 

2/2 

8,000 

9,300 

148 

0.0112 

273 

0.0082 

229 

10/10 

6,000 

9,300 

148 

0.0112 

273 

0.0082 

229 

10/10 

8,000 

9,300 

148 

0.0112 

273 

0.0082 

229 

10/10 

8,500 

9,300 

148 

0.0112 

273 

0.0082 

229 

10/10 

9,000 

9,300 

BH3 

55  ± 55 

150 

0.0136 

294 

0.0092 

246 

2/2 

5,000 

3,600 

150 

0.0136 

294 

0.0092 

246 

2/2 

6,000 

3,600 

BH4 

80±  30 

150 

0.0136 

294 

0.0048 

140 

2/2 

> 64,000 

16,700 

150 

0.0136 

294 

0.0048 

140 

2/2 

> 64,000 

16,700 

BH5 

70±  30 

148 

0.0112 

273 

0.0048 

140 

2/2 

185,000 

32,900 

148 

0.0112 

273 

0.0048 

140 

2/2 

> 185,000 

32,900 

BH9 

72.5±37.5 

150 

0.0136 

294 

0.0061 

175 

2/2 

72,400 

8,700 

150 

0.0136 

294 

0.0061 

175 

2/2 

96,000 

8,700 

BH7 

50±  50 

148 

0.0112 

273 

0.0082 

229 

105/2 

30,000 

36,900 

148 

0.0112 

273 

0.0082 

229 

105/2 

55,000 

36,900 

BH6 

55  ± 55 

150 

0.0136 

294 

0.0092 

246 

105/2 

8,000 

16,300 

150 

0.0136 

294 

0.0092 

246 

105/2 

9,000 

16,300 
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(13) 


log  N = -6.0796  - 4.8057  log  Ae  - 0.013647  (a^  + aie)  o (Ae)) 


provides  a method  for  evaluating  o^.  Equation  (13)  was  used  to  find  a value  of  for  each 
of  the  major  cycle,  rapid  (i.e.,  no  dwell  time)  test  conditions  (conditions  BHl,  BH2,  and 
BH3).  The  average  value  was  found  to  be  —160  ksi.  This  value  is  reasonably  close  to  the 
compressive  strength  of  the  material,  and  consequently  seems  to  be  a reasonable  value  for 
a residual  stress.  The  dwell  test  data  (conditions  BH6  and  BH7)  were  not  included  for  cal- 
culating Oq  in  order  to  preclude  the  possibility  of  relaxation  modeling  inaccuracies  affecting 
the  Oq  determination.  The  subcycle  conditions  BH4,  BH5,  and  BH9  were  not  included  for 
reasons  that  will  be  discussed  shortly. 

With  Oq  established,  the  relaxation  model  discussed  in  Appendix  C was  used  to  calculate  the 
mean  stresses  during  the  dwell  tests.  For  the  simple  cycle  tests,  the  relaxation  of  the  mean 
stress  can  be  written  as 

OR  = |oo+a(a)--i  a(Ae)^l  + ^14) 
and  the  instantaneous  mean  stress  can  be  written  as 

®m  = ®m  ^ a(Ae)^ ) (15) 

where  ~ ®o  — Y 

a = 0.002  hr"’ 


Figure  31  shows  the  relaxation  of  the  mean  stresses  for  the  simple  cycle  BH  test  conditions 
as  calculated  from  (15).  The  minimum  allowed  stress,  a^,  as  discussed  in  Appendix  C was 
taken  to  be  -160  ksi.  Substituting  (IS)  into  the  mean  stress  term  of  (13)  gives  a life  equation 
for  both  rapid  and  dwell  BH  test  conditions  as: 


log  N = -6.0796  - 4.8057  log  Ae 


-0.013647 


--.(Art  I 


t2  + ti) 


->] 


(16) 


Since  T is  a function  of  N,  an  iterative  procedure  was  used  to  find  the  predicted  life.  Result- 
ing predicted  f igue  lives  are  listed  in  Table  7,  and  the  resulting  correlation  is  shown  in 
Figure  32.  The  orrelation  of  both  rapid  and  dwell  major  cycle  data  is  good  with  the  single 
exception  of  one  specimen  tested  at  condition  BH I . This  specimen  failed  much  sooner 
than  expected  and  lies  outside  of  the  95  percent  confidence  band  established  by  the  SC 
data;  all  other  major  cycle  data  lies  well  within  the  confidence  band.  This  single  data  point 
is  suspected  as  a statistical  “outlier”.  Figure  33  compares  the  Weibull  distributions  resulting 
from  the  correlation  model  to  the  hypothetical  “perfect”  model  for  the  SC  and  the  major 
cycle  BH  data.  The  single  data  point  suspected  as  an  outlier  is  indicated  in  the  plot  but  was 
not  included  in  determining  the  slope  of  the  distribution.  Its  position  suggests  even  more 
strongly  that  it  is  an  outlier.  Once  again,  the  quality  of  the  correlation  model  is  judged  to  be 
good  for  the  major  cycle  data  based  on  the  results  shown  in  Figures  32  and  33. 


Figure  32  Correlation  of  Phase  II  SC  and  BH  Data 
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Figure  33  Phase  II  Combined  Data  Set  Weibull  Plots 


The  principal  influence  in  correlating  the  rapid  cvcle  test  points  is  the  recognition  of  the  sur- 
face layer  since  very  little  relaxation  takes  place  for  those  data  points  (see  Figure  31).  It  is 
clear  that  a surface  layer  does  exist  in  the  bolt-hole  specimens  and  that  its  behavior  is  mod- 
eled accurately  by  the  same  algorithm  as  employed  for  the  PWA  1216  bolt-hole  specimens. 
The  correlation  of  the  dwell  data  is  due  to  the  recognition  of  both  surface  layer  stresses  and 
their  subsequent  relaxation.  It  was  concluded  that  the  relaxation  model  incorporated  the 
principal  features  of  the  relaxation  process.  This  conclusion  was  further  supported  with  the 
cumulative  damage  testing. 

Figure  32  shows  that  the  subcycle  BH  test  conditions  are  consistently  underpredicted.  That 
is,  the  correlation  model  assigns  more  fatigue  damage  to  these  conditions  than  is  actually 
experienced.  This  was  definitely  not  the  case  for  the  PWA  1216  data.  This  comparison  sug- 
gests that  a material  dependent  mechanism  is  active  and  is  associated  with  subcycle  activity. 
Further,  comparison  of  local  stress/strain  conditions  for  the  major  cycle  and  subcycle  test 
conditions  (Figure  29)  indicates  that  mechanism  is  probably  associated  with  the  strain 
range.  This  behavior  is  suggestive  of  a fatigue  threshold  mechanism.  Regardless  of  the  mech- 
anism involved,  cumulative  damage  testing  using  these  subcycles  clearly  showed  that  sub- 
cycles contributed  to  damage  accumulation  when  used  in  combination  with  a major  cycle. 
Furthermore,  subcycles  can  never  exist  as  isolated  events  in  the  usage  of  a real  component. 
The  double-damage  cumulative  damage  model  discussed  in  Section  VI  fully  accounts  for  the 
subcycle  damage. 
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SECTION  VI 

CUMULATIVE  DAMAGE  LIFE  MODEL 
6. 1 THE  DOUBLE-DAMAGE  MODEL 

A cumulative  fatigue  damage  model  was  developed  using  the  double-damage  effect  described 
by  Manson^ . This  model  is  most  easily  illustrated  for  a block-loading  test  such  as  shown 
in  Figure  14.  For  such  a test,  the  familiar  linear  Miner’s  rule  assesses  the  damage  done  during 
each  block  of  loading  to  be 


Di  = ni/Ni  (16) 

where  Dj  is  the  relative  damage  due  to  ni  cycles  at  a load  condition  which  has  a simple 
cycle  life  of  N j . 

The  double-damage  effect  recognizes  that  the  entire  life  of  the  part  is  composed  of  a crack 
initiation  phase  and  a crack  propagation  phase.  For  a given  value  of  Dj,  the  crack  size  (follow- 
ing initiation)  developed  during  high  stress  amplitude  loading  such  as  the  a - loading  in  Figure 
14  is  larger  than  the  crack  that  would  have  been  developed  during  low  stress  amplitude  load- 
ing. Thus,  a low  amplitude  loading  block  following  a high  amplitude  loading  block,  during 
which  initiation  has  occurred,  will  result  in  failure  sooner  than  would  be  predicted  by  the 
linear  Miner’s  rule  damage  model 

Linear  damage  modeling  predicts  (for  the  block  loading  test  shown  in  Figure  14)  that  the 
cycles  to  failure  in  the  ^-loading  block,  np,  following  na  cycles  at  the  a - loading  block  is 
found  from  the  Miner’s  rule  equation 


na 

Na  N^ 


(17) 


If,  however,  a crack,  a,  has  been  initiated  at  the  end  of  the  a - loadmg  block,  then  the  double- 
damage model  implies  that  np  is  governed  by  fracture  mechanics. 

The  crack  growth  law  is  needed  to  derive  an  expression  for  np*,  the  cyclic  life  remaining 
under  p - loading  as  predicted  by  fracture  mechanics.  The  crack  growth  law  is  given  by 

~ = C AK"  (18) 
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where 


da 

— — = crack  growth  rate 
dN 

AK  = alternating  stress  intensity  factor 
n = crack  growth  exponent 
C = crack  growth  constant 


The  alternating  stress  intensity  factor  AK  must  include  the  effects  of  local  stress  ratio,  R 
(-  Omin/®max)-  For  the  materials  considered  in  this  contract,  the  R-ratio  model  can  be  in- 
corporated in  the  equation  for  AK  as 


where 


a = crack  size 

Y = crack  shape  function 

Ao  = local  stress  range  [=  a(Ae)] 

b,  d = material  constants  derived  from  R-ratio  crack 
growth  tests 


Integration  of  the  crack  growth  law  (18)  gives  an  expression  for  tifi*,  the  number  of  cycles 
required  to  grow  the  crack,  a to  a final  length  of  1/32  inch  length  under  |3  - loading. 


* s 


C k/T  Y 


d-R^ 


1 n 

P' 


Aom  (n/2  - 1) 


(20) 


A 

If,  instead,  the  crack  a had  been  propagated  under  a - loading  conditions,  the  remaining  life 
would  simply  have  been  equal  to  - n^.  And  integration  of  (18)  would  now  result  in  the 
expression 


Na-na 


1 


-n/2 


Combining  (20)  and  (21), 


^No(  - no! 


) 


(21) 


(22) 
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Manipulating  (22),  the  following  modifii*  ‘ Miner’s  rule  is  obtained: 


Thus,  double-damge  effects  have  the  influence  of  modifying  the  Miner’s  linear  damage 
rule.  It  can  be  shown  that  the  double-damage  debit  factor,  Da/j,  is  less  than  one  for  low 
amplitude  loading  following  high  amplitude  loading  (a  “High-Low”  test)  and  is  greatei 
than  one  for  high  amplitude  loading  following  low  amplitude  loading  (a  “Low-High” 
test).  For  (23)  to  apply,  it  is  only  necessary  that  a coherent  crack  be  initiated  at  some  time 
during  a - loading. 


While  the  double-damage  model  can  be  written  in  the  form  of  a modified  Miner  equation, 
(23),  it  should  be  remembered  that  it  is  actually  a fracture  mechanics  approach  which 
recognizes  the  exhaustion  of  the  initiation  life  during  the  first  loading  block.  This  fact  is 
further  illustrated  by  writing  (23)  in  an  expanded  form  and  noting  that  the  simple  cycle 
life  of  the  second  loading  block,  Njj,  does  not  enter  into  the  calculation  of  np*  . 


(24) 


Stated  yet  another  way,  the  double-damage  equation  (24)  simply  evaluates  the  crack  growth 
during  the  second  loading  by  evaluating  the  severity  of  crack  growth  conditions  in  the  second 
loading  block  relative  to  the  loading  which  exhausted  the  initiation  life. 

After  a crack  has  been  initiated,  it  propagates  both  along  the  surface  of  the  bolt  hole  as  well 
as  into  the  subsurface  region.  It  is  reasonable  to  assume  that  the  propagation  of  the  crack 
is  governed  essentially  by  the  more  positive  R-ratio  subsurface  region.  For  this  reason,  sub- 
surface stress  ranges  and  R-ratios  are  used  in  the  double-damage  equation  (24). 

As  stated  above,  for  (24)  to  apply,  it  is  only  necessary  that  a crack,  a,  exist  at  the  end  of 
the  a - loading  block  in  a block  - loading  type  of  test.  If  a does  not  exist  at  that  time. 

Miner’s  rule  initiation  life  exhaustion  damage  (16)  continues  to  accumulate.  In  such  a case, 
for  a two-block  block  loading  test,  the  double-damage  rule  would  reduce  to  the  simple 
Miner’s  Rule. 

For  more  complicated  mission  type  loading  such  as  the  sequenced  test  shown  in  Figure 
14,  it  is  clear  that  the  initiation  life  is  not  exhausted  by  just  one  of  the  various  type  of  cycles 
but  is  instead  exhausted  by  major  and  minor  cycles  acting  together.  As  such,  a does  not 
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exist  until  after  some  unknown  number  of  missions.  As  stated  previously,  ^rior  to  the  exist- 
ence of  a,  the  initiation  life  exhaustion  is  given  by  (16)  and  not  until  after  a exists  is  it 
appropriate  to  apply  the  double-damage  rule.  The  assumption  is  made  in  this  program 
that  all  minor  cycles  contribute  to  the  propagation  life  exhaustion.  Thus,  some  flight-by- 
flight subcycles  experienced  early  in  the  component  life  will  not  see  the  small  crack  initiated 
by  the  major  cycles  and  the  assumption  will  be  conservative.  For  very  long  lived  compon- 
ents, this  assumption  becomes  increasingly  more  conservative. 

6.2  PHASE  1 - CUMULATIVE  DAMAGE  TESTS 
6.2.1  TFR  Specimen  Testing 

One  gang  of  four  bolt-hole  specimens  was  used  to  determine  the  fatigue  damage  done  during 
the  terrain  following  radar  (TFR)  portion  of  an  actual  mission.  The  specimen  nominal  stress 
spectrum  shown  in  Figure  34  was  developed  by  considering  a random  sequence  of  ±7  percent 
(maximum)  PLA  excursions  for  the  FlOO  second-stage  fan  bolt-hole  location.  The  loading 
was  repeated  continuously  until  390,000  subcycles  were  imposed,  corresponding  to  roughly 
6,500  fleet-typical  flights  of  the  TF30  engine  in  the  F-1 1 1 . Testing  was  terminated  at  this 
point  and  subsequent  testing  of  one  of  the  TFR  specimens  at  test  condition  BH 1 (50  ± 50  ksi 
nominal  stress)  showed  that  the  1/32  inch  crack  life  was  no  different  from  virgin  specimens 
tested  under  the  same  loading  (see  Figure  35).  It  was  therefore  concluded  that  the  TFR 
specimens  could  be  considered  to  be  virgin  specimens.  One  of  these  specimens  was  subsequent- 
ly tested  at  simple  cycle  condition  BH4.  The  results  of  that  test  were  reported  in  SectionV. 

The  remaining  two  specimens  were  subsequently  tested  at  nominal  stresses  of  33  ± 33  ksi  for 
1 55,000  cycles  and  then  at  40  ± 40  ksi  until  failure.  A discussion  of  this  test  (labeled  BH54) 
is  included  in  the  next  section  (6.2.2). 


Figure  34  Terrain  Following  Radar  (ThR)  Simulation 
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LIFE  TO  1/32  INCH  CRACK  LENGTH  (CYCLES) 

Figure  35  Bolt  Hole  Specimens  Tested  at  50  ±50  Ksi  Nominal  Stress 
6.2.2  Cumulative  Damage  Specimen  Tests 


Phase  I BH  and  SC  cumulative  damage  test  conditions  and  resulting  lives  are  shown  in  Table 
8,  along  with  corresponding  Miner’s  Rule  (17)  and  double-damage  (24)  life  predictions. 

The  double-damage  debit  factor,  Da^,  used  in  these  predictions  is 


Do/J 


/ 1.64 -RA  4-84 

\1.64-Rc,/  \^0  ) 


(25) 


The  PWA  1216  room  temperature  crack  propagation  rate  exponent  (n  = 4.84)  and  the 
R-ratio  model  were  established  from  independent  testing. 


Tests  BH12  and  BH23  can  be  considered  “High-Low”  tests  (i.e.,  high  amplitude  loading 
preceeding  low  amplitude  loading)  and,  as  noted  previously,  the  double-damage  debit 
factors  for  these  tests  are  less  than  one.  Conversely,  tests  BH21  and  BH  32  are  “Low-High” 
tests  and  have  Dap's  greater  than  one.  Predicted  simple  cycle  lives  (as  discussed  in  Section 
V)  were  used  in  both  the  Miner’s  Rule  and  Double-Damage  calculations  shown  in  Table 
8.  Consequently  the  predicted  cumulative  damage  lives  are  influenced  by  the  simple  cycle 
life  model’s  ability  to  predict,  thus  making  an  evaluation  of  the  two  cumulative  damage 
models  more  difficult  .So,  in  order  to  allow  a more  accurate  evaluation,  the  average  simple 
cycle  lives  for  each  of  the  loading  blocks  were  first  used  in  the  Miner’s  Rule  and  double- 
damage equations.  The  Weibull  plots  of  these  results  are  shown  in  Figures  36a  and  36b. 

In  Figure  36a  it  is  clear  that  the  actual  lives  of  “High-Low”  tests  are  consistently  over-pre- 
dicted and  that  the  “Low-High”  tests  are  consistently  under  predicted.  However,  the  double- 
damage method  predicts  both  “High-Low”  and  “Low-High”  tests  uniformly  as  shown  in 
Figure  36b.  Thus,  double-damage  effects  are  clearly  displayed  and  are  predicted  by  the  double- 
damage equation  (24).  Figure  36c  shows  the  effect  of  using  predicted  simple  cycle  lives  in 
(24);  again,  the  tests  appear  to  be  uniformly  predicted. 
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Figure  36  Nonlinear  Cumulative  Damage  Effects 

A 

The  clear  existence  of  the  double-damage  effect  implies  that  a crack,  a,  existed  at  the  end 
of  the  first  loading  block  in  these  tests.  In  each  of  the  tests,  the  number  of  cycles  applied 
in  the  first  loading  block  was  approximately  50  percent  of  the  number  of  cycles  needed  to 
initiate  a small  (approximately  0.008  inch)  crack  at  that  loading  condition  as  established 
from  simple  cycle  testing.  Thus,  it  seems  reasonable  that  a small  sub-detectable  crack,  a, 
could  indeed  have  existed  at  the  end  of  the  first  loading  block. 

The  test  labeled  BH54  is  included  in  Table  8 for  completeness  even  though  it  is  judged 
not  to  be  directly  applicable  for  the  purposes  of  this  contract.  As  shown  in  Table  8, 

BH54  is  correlated  by  neither  the  Miner’s  Rule  nor  the  double-damage  equation,  but 
examination  of  the  calculations  shown  in  the  table  indicate  some  obvious  limits  in  the  pre- 
diction systems.  First,  the  calculated  simple  cycle  life  (3,400,000  cycles)  for  the  first 
loading  block  is  obviously  much  too  optimistic  since  the  specimens  failed  shortly  after  the 
first  155,000  cycles  had  been  applied.  This  over  prediction  is  not  surprising  since  the  data 
base  is  limited  to  much  lower  lives.  Secondly,  the  double-damage  calculation  predicts  that 
the  effect  of  is  to  enhance  the  life  of  the  second  loading  block.  That  is,  the  prior  load- 
ing is  predicted  to  have  caused  n^*  to  be  larger  than  the  simple  cycle  life  of  the  second 
block,  N/j.  Clearly,  this  implies  a limit  on  Da/3  for  “Low-High”  tests.  This  limit  has  been 
exceeded  somewhat  in  condition  BH32  as  well,  although  it  has  not  been  exceeded  as  sub- 
stantially as  exhibited  in  the  calculation  for  BH54.  Lastly,  such  a limitation  of  the  double- 
damage calculation  occurs  only  in  the  “Low-High”  type  of  loading  (where  Da/3>  ' 0).  This 
type  of  loading  never  occurs  in  component  usage.  For  this  reason,  the  implied  limit  in  the 
double-damage  calculation  does  not  have  to  be  considered  further  for  the  cumulative  dam- 
age predictions  considered  in  this  program. 
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The  correlation  of  the  cumulative  damage  specimen  tests  (including  the  single  SC  test) 
using  the  double-damage  approach  (25)  is  shown  in  Figure  37.  Values  of  np*  are  plotted 
for  the  block  loading  tests. 
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Figure  37  Correlation  of  Phase  I Simple  Cycle  and  Cumulative  Damage  Specimen  Testing 


6.2.3  Ferris  Wheel  Tests 


Verification  of  the  fatigue  life  prediction  models  described  in  previous  sections  was  achieved 
by  tests  of  four  fan  test  disks.  The  disk,  shown  in  Figures  38  and  39  was  designed  to  be  LCF 
life  limited  in  the  bolt  hole.  The  elastic  stress  concentration  factor  (Kt  = 2.427)  was  deter- 
mined from  extensive  finite  element  and  finite  difference  analyses  and  verified  with  strain 
gage  surveys.  Disk  testing  was  performed  in  a “Ferris  Wheel”  (F/W)  which  simulates  engine 
operating  stresses  through  radial  loading  applied  at  the  disk  rim. 

The  generic  missions  established  from  the  mission  surveys  reported  in  Section  II  were 
used  as  the  basis  for  the  disk  testing.  Actual  disk  nominal  stress  levels  and  the  resulting  lives 
are  shown  in  Table  9.  The  major  cycle  local  stresses  and  strains  of  each  Ferris  Wheel  disk 
test  are  shown  in  relation  to  the  BH  specimen  test  condition  in  Figure  40.  Three  of  the  disk 
tests  were  run  to  final  failure  but  F/W4  test  was  suspended  with  no  bolt  holes  cracked.  All 
three  disks  tested  to  completion  failed  from  a single  bolt  hole;  postfailure  examinations 
revealed  no  crack  indications  in  any  of  the  remaining  holes. 
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PHASE  I FERRIS  WHEEL  TEST  CONDITIONS  AND  RESULTS 


missions  C< 


-8 


Figure  40  Specimen  and  Ferris  Wheel  Surface  Layer  and  Subsurface  Response 

The  life  prediction  methods  described  in  previous  sections  have  been  developed  for  pre- 
dicting typical  or  mean  LCF  life  behavior,  but  because  each  test  disk  could  have  failed 
from  any  one  of  twenty  bolt  holes,  life  predictions  for  the  test  disks  were  made  for  the  appro- 
priate probability  of  failure  rate.  Based  on  the  probability  distribution  of  actual  versus 
average  lives  for  the  BH  specimens,  the  lowest  lived  of  twenty  holes  would  fail  at  65  percent 
of  the  mean  predicted  LCF  life. 

The  disk  bolt  holes  were  manufactured  using  machining  parameters  which  duplicated  those 
used  in  the  BH  specimen  machining.  However,  the  bolt  holes  in  the  second  and  third  disks 
were  found  to  have  machining  marks  atypical  of  the  BH  specimens  and  were  subsequently 
reoperated  by  honing  to  a maximum  diametral  increase  of  0.008  inch.  Thus,  the  surface 
layer  of  the  second  and  third  disks  were  not  the  same  as  for  the  first  and  fourth  disks. 

The  double-damage  model  discussed  in  previous  sections  was  used  to  predict  the  number  of 
missions  required  to  fail  the  first  hole  in  each  disk.  Due  to  the  very  low  amplitude  of  loading 
in  Mission  A of  F/W3,  this  mission  was  assumed  to  contribute  only  to  the  exhaustion  of 
initiation  life;  no  double-damage  effects  are  assumed  to  be  present  in  Mission  A.  Figure  41 
shows  how  the  correlation  of  actual  and  predicted  missions  varies  with  the  assumed  values 
of  surface  layer  machining  stress,  Oq.  It  is  clear  that  F/W2  and  F/W3  behave  differently  than 
F/Wi . Furthermore,  the  value  of  o©  (=  -133  ksi)  which  correlates  F/Wj  is  nearly  identical 
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to  oq  determined  for  the  BH  specimens.  This  correlation  is  expected  since  these  parts  were 
machined  in  the  same  way.  However,  the  value  of  a©  (=  —71  ksi)  which  best  correlates  F/W2 
and  F/W3  reflects  the  fact  that  part  of  the  surface  layer  was  removed  by  subsequent  honing. 
This  is  viewed  as  clear  evidence  of  the  existence  of  the  surface  layer  and  firmly  establishes 
validity  of  the  surface  layer  algorithm.  The  predicted  lives  for  each  of  the  F/W  tests  is  shown 
in  Table  9. 


Figure  41  Surface  Residual  Stress  for  the  Second  and  Third  Disks  Shows  the  Effect  of 
Partial  Surface  Layer  Removal 


Overall  correlation  of  all  specimen  and  Ferris  Wheel  test  data  is  shown  in  Figures  42  and 
43.  This  correlation  is  judged  to  be  extremely  good  considering  the  wide  variety  of  speci- 
mens, components,  and  loading  conditions  considered. 


Figure  42  Correlation  of  Phase  I Specimen  and  Ferris  Wheel  Data 
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Figure  43  Phase  I Overall  Correlation 


6.3  PHASE  II  - CUMULATIVE  DAMAGE  TESTS 
6.3.1  Cumulative  Damage  Specimen  Tests 

The  PWA  1057  cumulative  damage  investigation  involved  several  block  loading  tests  and 
one  sequenced  test.  The  actual  test  conditions  and  resulting  lives  are  shown  in  Table  lO. 

The  Miner’s  Rule  and  double-damage  life  calculations  for  each  test  condition  are  also  shown 
in  the  table.  The  double-damage  debit  factor  used  in  these  predictions  is 

/I  077  -Rg\  2 

The  PWA  1057  WO^F  crack  propagation  rate  exponent  (n  = 2.0)  and  the  R-ratio  model 
were  established  from  independent  testing. 

Consider  first  those  tests  involving  only  rapid  cycle  testing;  BH32-2,  BH29-1 , and  BH29-2. 

In  test  BH32-2  the  Miner’s  Rule  life  calculation  is  the  better  of  the  two  methods  for  correl- 
ating the  life.  This  is  reasonable  since  it  is  relatively  certain  that  a crack,  a,  was  not  initiated 
during  the  first  loading  block  (a  - loading).  It  required  approximately  2000  cycles  at  a - load- 
ing level  to  initiate  a “pin-point”  crack  (approximately  0.008  inch)  in  the  simple  cycle  test- 
ing. As  stated  in  Section  6.1,  the  double-damage  concept  reduces  to  Miner’s  Rule  for  such 
cases.  Tests  BH29-1  and  29-2  illustrate  again  that  the  double-damage  model  is  appropriate. 
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While  the  Miner’s  Rule  calculation  appears  at  first  glance  to  be  comparable  to  the  double- 
damage calculation,  it  should  be  remembered  that  the  simple  cycle  life  prediction  for  np*  can 
too  low  by  nearly  an  order  of  magnitude  as  reported  in  Section  5.0.  Thus,  a Miner’s  Rule 
prediction  for  n^  using  an  accurate  simple  cycle  prediction  for  would  be  much  higher  than 
the  observed  np.  But,  as  discussed  in  Section  6.1,  the  double-damage  equation  for  n^*  can 
be  written  in  the  form  (22)  in  which  does  not  appear.  So  the  double-damage  calculation 
leads  to  a more  accurate  answer. 

It  is  clear  then  that  while  the  subcycle  loading  (P  - loading)  produces  essentially  no  damage 
in  itself,  it  does  contribute  significantly  to  propagation  life  exhaustion. 

Test  condition  BH27  demonstrates  that  once  the  initiation  life  has  been  exhausted  by  prior 
loading,  the  dwell  cycle  contributes  to  propagation  life  exhaustion  in  the  manner  described 
by  the  double-damage  equation.  It  was  therefore  inferred  that  the  principal  influence  of 
the  dwell  cycle  is  to  enhance  the  initiation  life. 

Subsequent  block  loading  tests  BH62,  BH73-1,  and  BH73-2  were  conducted  to  investigate 
the  influence  of  dwell  cycles  on  initiation  life.  This  was  done  by  selecting  the  number  of 
cycles  in  the  first  loading  block,  n^,  such  that  a crack  did  not  exist  at  the  end  of  the  a - load- 
ing block.  Thus  for  these  tests,  the  Miner’s  Rule  life  prediction  is  appropriate.  The  principal 
effect  of  prior  dwell  loading  on  subsequent  loading  was  found  to  be  a reduction  of  the  mean 
stress  levels  over  which  subsequent  cyclic  activity  occurred.  This  stress  condition  is  reflected 
in  the  simple  cycle  lives,  Np,  for  the  0 - loading  condition.  For  the  double-damage  calculation, 
the  relaxation  of  both  surface  layer  (affecting  Np  values)  and  the  subsurface  region  (affecting 
the  crack  growth  calculation)  was  recognized.  The  relaxation  of  the  mean  stress  of  the  sub- 
surface region  is  shown  in  Figure  44. 


Figure  44  Subsurface  Mean  Stress  Relaxation 
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In  a similar  manner,  for  the  sequenced  test,  BH73S,  the  effect  of  the  dwell  period  was 
taken  to  be  the  relaxation  of  the  mean  stress  level  of  the  rapid  cycles.  In  the  double-damage 
calculation,  the  propagation  behavior  of  the  dwell  subcycle  was  debited  relative  to  the  prin- 
cipal life  exhausting  cycle  (the  50  ± 50  ksi  rapid  cycle).  The  relaxation,  both  surface  and  sub- 
surface occurs  over  a period  of  several  hours  and  eventually  reaches  a fully  relaxed  condition 
as  discussed  in  Appendix  C.  The  simple  cycle  lives  and  double-damage  debit.  Dap  , shown 
in  Table  10  were  calculated  for  the  fully  relaxed  condition,  but  it  was  calculated  that  nearly 
10000  missions  would  be  required  to  achieve  the  fully  relaxed  condition.  It  is  apparent  that 
if  a partially  relaxed  condition  were  to  be  used  as  the  basis  for  calculation,  lower  lives  would 
be  predicted.  This  result  suggests  an  iterative  process  for  determining  the  Nj  and  Dij  for 
sequenced  tests.  From  Table  10,  it  appears,  once  again,  that  the  double-damage  method 
provides  the  better  correlation. 

Figure  45  shows  the  overall  correlation  of  the  cumulative  damage  tests.  Total  lives  (n^^  + na) 
are  plotted  here.  The  correlation  is  well  within  the  95  percent  confidence  band  established 
by  the  SC  data  alone. 

6.3.2  Ferris  Wheel  Tests 

Verification  tests  were  performed  on  Ferris  Wheel  disks  at  a constant  temperature  of 
900‘’F.  Figures  46  and  47  show  the  actual  disks  mounted  in  the  test  rig  with  the  heating 
mechanisms  in  place.  Figure  48  shows  the  dimensions  of  the  test  disk.  This  shape  contains 
the  major  features  of  a typical  turbine  disk.  The  elastic  stress  concentrations  factor  for  the 
disk  was  determined  to  be  2.33. 

Actual  test  loading,  based  on  the  generic  missions  of  Section  II,  is  shown  in  Table  1 1.  Each 
disk  was  subjected  to  a room  temperature  calibration  run  prior  to  testing  to  ensure  that 
proper  loading  conditions  were  being  achieved.  For  tests  F/W  j and  F/W^,  the  room  tempera- 
ture calibrations  were  made  using  a maximum  nominal  stress  equal  to  that  experienced  dur- 
ing 900°F  testing.  However,  in  the  room  temperature  calibration  of  F/W2,  the  maximum 
stress  was  16  percent  higher  than  the  maximum  stress  experienced  during  subsequent  cycling. 
The  effects  of  these  room  temperature  cycles  were  accounted  for  in  the  prediction  of  the 
F/W  test  results.  Figure  49  shews  schematically  how  the  effect  of  the  room  temperature 
loading  was  calculated.  Simply  stated,  the  room  temperature  monotonic  and  hysteresis 
curves  are  used  to  calculate  an  initial  stress  and  strain  state  (point  B in  Figure  49)  for  sub- 
sequent high  temperature  loading.  The  high  temperature  loading  continues  from  point  B 
along  the  ^OO^F  hysteresis  curve  for  both  the  surface  layer  and  surface  region. 

The  surface  layer  algorithm,  using  Oq  =-160  ksi  as  determined  from  BH  tests,  and  the  relaxa- 
tion algorithm  (Appendix  C)  were  used  to  calculate  the  local  bolt-hole  stresses  and  strains. 
The  double-damage  equation  was  used  in  the  life  prediction.  None  of  the  tests  reached  a 
fully  relaxed  condition  prior  to  the  end  cf  testing. 

In  all  cases  the  simple  cycle  lives,  Nj,  and  the  double-damage  debits,  Dij,  for  a given  mission 
were  calculated  using  the  stresses  occuring  at  the  end  of  that  mission.  For  the  Final  mission 
in  a test,  the  relaxation  was  calculated  at  the  number  of  missions  corresponding  to  first 
cracking.  All  life  predictions  were  made  for  typical  bolt-hole  behavior,  i.e.,  half  of  the  holes 
in  the  disk  showing  1/32  inch  cracks.  In  all  three  tests,  a sufficient  portion  of  the  bolt  holes 
were  cracked  to  establish  the  actual  typical  lives.  Typical  lives  are  shown  in  Table  1 1. 
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PHASE  11  CUMULATIVE  DAMAGE  SPECIMEN  TEST  CONDITIONS  AND  RESULTS 


Figure  4 6 Phase  H Perris  Wheel  Rig  ( 77-44 1-4 008- A ) 


Figure  45  Correlation  of  Phase  II  Simple  Cycle  and  Cumulative  Damage  Specimen  Testing 
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Figure  47  Phase  U Ferris  Wheel  Disk  in  Test  Rig  ( 77-44 J-4008-B) 


Figure  48  Dimensions  of  Phase  II  Ferris  Wheel  Disk 
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Figure  49  Local  a,  e Response  to  Room  Temperature  Loading 

In  disk  tests  F/W2  and  F/W^,  the  loading  was  changed  from  one  type  of  mission  to  another 
relatively  early  in  the  life  of  the  disk,  thus  creating  the  possibility  that  the  initiation  life  was 
not  exhausted  prior  to  the  change  in  mission  type.  As  stated  previously,  the  double-damage 
debit,  Dij,  is  properly  applied  only  to  those  cycles  occuring  after  the  initiation  life  has  been 
exhausted.  In  the  case  of  disk  test  F/W^,  the  assumption  that  the  initiation  life  was  exhausted 
(i.e.,  a existed)  during  the  first  applications  of  Mission  C leads  to  a double-damage  predic- 
tion that  the  disk  would  have  failed  during  Mission  B loading.  This  assumption  obviously 
was  not  the  case.  On  the  other  hand,  if  the  1300  Mission  C and  the  2500  Missions  B contrib- 
uted only  to  initiation  damage,  the  double-damage  model  would  predict  an  additional  3250 
Missions  C until  failure  of  half  of  the  bolt  holes  which  agrees  very  well  with  actual  life. 

It  is  clear  that  if  1 300  Missions  C were  insufficient  to  exhaust  the  initiation  life  in  F/W^,  then 
certainly  only  850  Missions  A would  also  not  exhaust  the  initiation  life  in  F/W-,.  The  double- 
damage life  model  then  predicts  that  an  additional  4430  Missions  C would  be  needed  to  fail 
50  percent  of  the  bolt  holes.  This  prediction  also  agrees  very  well  with  the  actual  typical  life 
as  shown  in  Table  1 1. 

Since  only  one  type  of  mission  was  used  in  disk  test  F/W  ].  all  subcycles  were  assumed  to 
contribute  to  the  propagation  of  the  crack,  a.  ('onsequcntly.  the  predicted  life  is  conservative. 
But  again,  the  agreement  between  actual  and  predicted  lives  is  good. 
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Figure  50  shows  the  overall  correlation  of  all  Phase  II  specimen  and  Ferris  Wheel  tests.  The 
• correlation  of  all  the  cumulative  damage  tests  is  within  the  95  percent  confidence  band  de- 

termine from  the  SC  data  alone.  The  Weibull  distribution  of  the  data  (subcycle  data  not  in- 
cluded as  discussed  previously)  is  shown  in  Figure  5 1 . This  overall  correlation  is  judged  to 
be  very  good  indicating  that  the  various  components  of  the  prediction  model  accurately  de- 
scribe the  damage  processes  that  occurred. 

i 6.4  CUMULATIVE  DAMAGE  PREDICTION  COMPUTER  PROGRAM 

! A computer  program  was  written  for  the  prediction  of  disk  cumulative  damage  life  exhaus- 

, tion  using  the  methods  discussed  in  this  report.  The  program  recognizes  the  existence  of  a 

j surface  layer  and  the  relaxation  of  stresses  caused  by  dwell  periods.  The  cumulative  damage 

effects  are  modeled  by  the  double-damage  method  described  above.  The  assumption  is  made 
that  all  minor  cycles  contribute  to  propagation  life  exhaustion.  The  computer  program  is  | 

discussed  further  in  Appendix  E.  Computer  program  documentation  is  contained  in  Ref-  * 

erences  1 and  2.  j 


% failed 


Figure  50  Correlation  of  Phase  II  Specimen  and  Ferris  Wheel  Data 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 
7.1  LCF  LIFE  EXHAUSTION  MODELS 

Prediction  of  low  cycle  fatigue  (LCF)  damage  due  to  complex  load  history  is  an  active  area 
of  research^’®  and  is  not  limited  to  aerospace  applications’.  A major  conclusion  of  the  cur- 
rent work  is  that  most  of  the  published  fatigue  results  (such  as  those  presented  in  Dowling 
et  al*®)  are  inappropriate  for  the  current  research  effort  due  to  the  use  of  balanced  strain 
testing  which  neglects  mean  stress  effects.  Further,  it  is  concluded  that  no  acceptable  design 
model  can  be  generated  that  does  not  reflect  differences  due  to  specimen  type  and  machin- 
ing effects. 

The  approach  used  in  the  reported  research  and  development  program  was  the  local  strain 
approach,  as  discussed  by  Dowling  et  al;  that  is,  unnotched,  strain  controlled  specimens 
are  tested  at  conditions  of  strain  range  and  mean  stress  that  represent  those  of  the  loaded 
notch.  The  notch  conditions  were  necessarily  modeled  in  a manner  that  accounted  as 
fully  as  possible  for  nonlinear  response,  in  a cost-effective  manner.  However,  it  was  clearly 
concluded  that  notched  specimens,  representative  of  the  actual  component  notch  including 
all  machining  operations,  must  be  tested  in  order  to  predict  component  life.  Further  conclu- 
sions on  machining  effects  are  discussed  in  Section  7.2. 

Simple  correlation  models  such  as  used  in  the  current  effort  can  successfully  fit  LCF  data 
taken  over  a wide  range  of  loading  and  specimen  conditions  if  proper  care  is  taken  to  define 
the  local  strain  and  mean  stress  conditions.  The  approach  used  emphasizes  automated  cor- 
relation methods,  as  opposed  to  engineer-interpreted  design  curves.  With  care,  some  of  the 
“scatter”  between  actual  lives  and  the  simple  correlation  model  can  be  reduced  by  local  ad- 
justments in  the  life  curves.  However,  the  risk  that  the  adjustment  could  be  in  error  for  other 
load  conditions  is  perhaps  as  great  as  the  original  error.  It  is  concluded  that  the  simple  correla- 
tion model  is  as  good  an  approach  for  design  as  hand  drawn  life  curves,  for  the  data  obtained 
in  the  current  program. 

The  life  prediction  model  is  based  on  correlating  the  number  of  LCF  cycles  to  initiate  a 1 /32 
inch  surface  length  crack  in  the  specimens  and  components.  It  was  concluded  that  while 
some  propagation  of  a coherent  crack  was  involved  in  the  LCF  initiation  life,  the  effects  of 
different  stress  (strain)  gradients  was  not  important.  Further,  it  was  concluded  that  all  lives 
(with  a possible  exception  for  very  short  LCF  life  data)  included  demonstrable  initiation 
and  propagation  behavior.  While  further  research  into  microcrack  growth  modeling  is 
recommended,  the  use  of  a fracture  mechanics  model  for  the  LCF  life  data  obtained  is  not 
recommended'* . 


7.2  SURFACE  LAYER  MODELS  AND  MACHINING  EFFECTS 


Clearly,  one  of  the  most  important  conclusions  from  the  current  research  and  development 
effort  was  the  need  to  systematically  account  for  a surface  layer  of  previously  worked  mate- 
rial in  the  tested  notches.  White  it  is  generally  accepted  that  both  residual  stress  and  surface 
cyclic  work  hardening  are  important  effects'^,  the  authors  are  not  aware  of  a simple  model 
for  these  effects,  comparable  to  the  one  successfully  adopted  in  the  current  effort. 

Scatter  in  LCF  life  data  is  a major  influence  on  minimum  design  life  calculations  for  a safe- 
life  design  methodology.  Control  of  scatter  to  increase  the  minimum  cyclic  life  is  a critical 
area  for  future  research.  The  simple  surface  layer  model  developed  in  the  current  effort  is 
based  only  on  bolt  holes,  machined  under  carefully  controlled  conditions.  Further  research  is 
required  to  establish  better  models  for  the  nature  of  the  surface  layer  work  hardening,  the 
level  of  the  imposed  residual  stresses,  and  the  effects  of  different  machining  operations.  While 
such  an  effort  is  a major  undertaking,  the  potential  pay-off  appears  to  be  quite  attractive  in 
terms  of  increased  performance  and  reduced  cost  for  critical  rotating  structure  components. 

7.3  DWELL  (CREEP)  EFFECTS 

The  presence  of  time  dependent  behavior  in  turbine  disks  was  clearly  established  in  the  cur- 
rent effort.  All  of  the  elevated  temperature  testing  emphasized  dwell  events  at  high  stress 
levels  which  generally  benefit  life.  However,  the  analytical  model  adopted  for  dwell  effects 
on  mean  cyclic  stress  includes  the  potential  for  life  debits  due  to  dwell  at  low  stress  levels. 
Further  testing  to  verify  the  potential  life  debit  is  recommended. 

It  was  generally  concluded  that  no  acceptable  material  behavior  analytical  model  exists  for 
predicting  the  time-dependent  creep  effect  which  can  account  for  cyclic  loading  at  locally 
high  stress  levels.  Thus,  cyclic  test  data  had  to  be  obtained  which  could  be  used  to  develop 
a simple  exponential  decay  model  for  cyclic  dwell  stresses.  The  model  used  was  based  on 
data  at  one  stress  level  and  one  temperature.  Further  work  is  required  to  develop  a calibrated 
model  of  dwell  stress  relaxation  for  different  stresses  and  temperatures. 

Finally,  the  current  effort  did  not  involve  creep-fatigue  interaction  effects.  Wlrile  such  problems 
also  involve  time-dependent  plasticity,  the  temperature  levels  of  the  current  effort  are  not 
judged  to  be  sufficient  to  contribute  creep  damage  to  the  material  microstructure.  Major 
new  research  efforts  are  required  to  model  creep-fatigue  damage  at  a level  of  confidence 
comparable  to  that  achieved  in  the  current  effort. 

7.4  CUMULATIVE  DAMAGE  MODEL 

An  early  conclusion  in  the  research  and  development  effort  was  that  the  nature  of  LCF 
cycles  in  military  gas  turbine  engine  disks  contained  one  critical  simplification.  Essentially 
all  important  cyclic  stress  excursions  at  the  disk  bolt-holes  shared  a common  maximum 
stress.  Thus,  nonlinear  overload  damage  retardation  modeling  was  not  required.  This  conclu- 
sion will  generally  apply  to  all  disk  notch  locations  with  the  exception  of  those  seeing  major 
transient,  thermally  driven  excursions. 
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The  oa'y  nonlinear  damage  effect  that  remained  was  named  by  Manson'^  as  the  “double- 
damage” effect.  Simply  stated,  a stress  (strain)  cycle  does  damage  at  a different  rate  on 
cracked  material  than  on  uncracked  material.  The  data  generated  in  the  current  program 
clearly  supported  this  hypothesis.  A fracture  mechanics  based  model  was  developed  which 
accounts  for  this  effect  by  debiting  the  cyclic  initiation  life  of  low  amplitude  cycles  to 
account  for  the  cracking  damage  done  by  the  major  cycles.  While  somewhat  conservative, 
the  adopted  algorithm  was  found  to  give  good  nonlinear  damage  estimates  for  military  mis- 
sion modeling. 

7.5  DESIGN  MODELS  FOR  LCF  LIFE  EXHAUSTION 

The  model  developed  by  the  current  effort  is  judged  to  be  suitable  for  design  purposes  under 
a few  stipulated  limitations.  The  current  effort  was  calibrated  only  for  disk  bolt-holes,  al- 
though the  procedures  are  felt  to  be  generally  applicable.  Single  heats  of  the  two  materials 
and  single  machining  operations  were  included  in  the  current  statistical  base.  Actual  design 
implementation  of  the  current  model  must  include  the  greater  statistical  scatter  associated 
with  heat-to-heat  and  processing  variables.  Application  of  the  current  model  to  hot  section 
components  is  limited  to  zero  creep-fatigue  interaction  and  to  mission  cycles  represented  by 
nearly  constant  values  of  peak  tensile  cyclic  stresses.  Also,  calibration  of  the  dwell  model 
used  is  limited,  as  discussed  in  Section  7.3,  to  a single  temperature  and  dwell  stress  level. 

Use  of  the  developed  LCF  life  exhaustion  model  for  engine  disk  design  is  possible  for  new 
alloy  systems  where  data  generation  can  be  obtained  along  lines  similar  to  those  used  in  the 
reported  program.  Application  of  the  developed  model  to  old  materials  is  not  recommended 
without  the  development  of  the  necessary  new  data  base.  In  both  cases,  it  should  be  clearly 
recognized  that  the  development  of  a data  base  for  LCF  life  exhaustion  modeling  for  design 
purposes  is  a major  undertaking. 
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APPENDIX  A 

MATERIAL  CHARACTERIZATION 

Basic  material  characterization  tests  were  conducted  for  a single  heat  code  of  each  of  two 
disk  alloys;  Ti  (6Al-2Sn-4Zr-6Mo),  (PWA  1216),  heat  code  CAAZ  and  WASPALOY® 
(PWA  1057),  heat  code  XNNZ.  Monotonic  and  cyclic  uniaxial  stress-strain  tests  were  per- 
formed for  both  materials  using  the  close-loop,  servo-controlled  MTS  testing  machine.  The 
specimen  used  is  the  strain  controlled  specimen  shown  in  Figure  A-1 . Gage  section  strains 
were  controlled  within  ±0.001  in/in;  loads  were  measured  with  an  accuracy  of  ±30  lb.  In 
addition  cyclic  relaxation  and  standard  creep  tests  were  performed  for  PWA  1057  using 
the  specimens  shown  in  Figures  A-1  and  A-2  respectively. 


Figure  A-1  Strain  Controlled  Constitutive  Test  Specimen 


0.620  0.2S1 


Figure  A-2  Creep  Specimen 


I 

I. 

' The  monotonic  stress-strain  curves  are  shown  in  Figures  A-3  and  A-4.  No  strain  rate  sensi- 

I tivity  was  observed  over  the  range  of  strain  rates  of  interest;  0.01  in/in/min  to  0. 10  in/in/ 

min.  Yield  behavior  of  the  PWA  1057  was  discontinuous  at  900‘’F.  Figure  A-5  is  a tracing 
I of  an  actual  load-deflection  plot  showing  this  discontinuous  behavior.  These  slip  bursts 

! were  accompanied  by  accoustical  emissions  and  were  observed  only  during  monotonic 

. straining  greater  than  approximately  0.007  in/in.  The  load  drop  during  these  slip  bursts 

were  within  the  variation  of  the  monotonic  strengths  from  repeated  tests  and  were  therefore 
i not  modeled  explicitly  by  the  monotonic  curve  adopted  for  use  in  this  program. 


6,  Monotonic  strain  (%  strain) 

Figure  A-3  PWA  1216  Room  Temperature  Monotonic  Stress-Strain  Curve 


6, Monotonic  strain  (%  strain) 
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Deflection,  (in.) 


Figure  A-5  Load-Deflection  Tracing  Showing  Discontinous  Yielding  of  PWA  1057  at  900° F 

Incremental  step  tests  were  used  to  characterize  the  cyclic  stress-strain  behavior  of  the  two 
materials.  Cycling  was  strain  controlled  at  a fixed  strain  rate  between  selected  strain  limits 
until  stabilized  hysteresis  loops  were  observed.  Maximum  strain  limits  were  incrementally 
increased  while  maintaining  the  same  minimum  strain.  Figure  A-6a  illustrates  the  technique. 

In  both  PWA  1216  and  PWA  1057,  stabilized  hysteresis  loops  were  obtained  within  a few 
cycles.  Neither  material  exhibited  cyclic  hardening  or  cyclic  softening.  The  cyclic  stress- 
strain  behavior  was  characterized  by  plotting  the  shapes  of  the  hysteresis  curves,  i.e.,  curves 
A,  B,  and  C for  each  test.  This  is  illustrated  in  Figure  A-6b.  For  PWA  1216,  the  shape  of 
each  hysteresis  loop  formed  a single  hysteresis  curve  (Figure  A-7),  while  in  PWA  1057,  each 
strain  range  had  its  own  characteristic  curve  shape  (Figure  A-8).  The  loci  of  tips  of  such 
curves  were  taken  to  represent  the  cyclic  behavior  of  PWA  1057.  The  resulting  curves  are 
shown  in  Figure  A-9. 

Standard  creep  tests  of  PWA  1057  at  900°F  were  conducted  using  the  specimen  shown  in 
Figure  A-2.  The  results  of  these  tests  are  shown  in  Figure  A- 10. 

Two  cyclic  relaxation  tests  were  also  conducted  using  the  specimen  shown  in  Figure  A-1 . 
During  these  tests,  the  specimen  was  cycled  between  fixed  strain  limits  with  a 125  second 
“dwell”  time  imposed  at  the  maximum  strain  of  each  cycle.  While  total  stress  range  was 
unchanged  during  cycling  (neither  hardening  nor  softening)  the  peak  stress,  and  conse- 
quently the  mean  stress  was  observed  to  change  during  cycling.  These  results  are  shown  in 
Figure  A-1 1.  The  stress  relaxation  was  found  to  be  approximated  by  equation  A-1. 

a=  155e'0-®®2t 

where  a = the  current  stress  (KSI)  (A-1) 

t = time  spent  at  peak  stress  (hr) 


(b) 

Figure  A-6  Incremental  Step  Tests  Used  to  Establish  Hysteresis  Curve 


Figure  A-7  PWA  1216  Hysteresis  Curve  Shape 


Figure  A-8  PWA  1057  Hysteresis  Curves  Show  Dependence  Upon  Strain  Range 


Figure  A-9  FWA  1057  Hysteresis  Curve  Shapes 
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Figure  A-10  Results  ofPWA  l057  Cyclic  Relaxation  Test 


t,  Total  time  at  max  strain  (hrs) 


Figure  A-II  Results  of  PWA  1057  Cyclic  Relaxation  Test 
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APPENDIX  B 

SURFACE  LAYER  ALGORITHM 

The  following  procedure  for  calculating  local  surface  mean  stress  recognizes  the  influences 
of  residual  surface  stresses  due  to  machining  and  cold  working.  The  influence  of  residual 
stresses  was  found  to  be  of  prime  importance  in  developing  a fatigue  life  prediction  model 
that  encompassed  both  Strain  Controlled  (SC)  and  Bolt  Hole  (BH)  specimens. 

The  bolt  hole  region  can  be  viewed  as  a thin  preworked  surface  layer  having  an  initial  stress 
state  (Oq,  e^j)  lying  atop  a subsurface  region  of  virgin  material  (Figure  B-1).  Because  the  re- 
sidual machining  stresses  do  not  extend  appreciably  below  the  bolt  hole  surface,  the  bulk 
response  of  the  region  is  not  significantly  altered  by  the  prestressed  surface  layer  and  can 
therefore  be  modeled  by  the  modified  Neuber  approach,  as  given  in  Section  4.0.  Compati- 
bility requires  that  both  surface  and  subsurface  strain  excursions  are  identical.  However, 
the  stress-strain  response  of  the  two  regions  will  be  different  since  the  surface  layer  has 
already  experienced  large  plastic  strains  during  machining,  and  is  therefore  responding  ac- 
cording to  the  hysteresis  curve  shape,  whereas  the  subsurface  region  is  not. 

Consequently,  during  the  loading  portion  of  the  first  cycle  (A  to  B in  Figure  B-1),  as  the 
subsurface  region  strains  to  an  amount  e (determined  by  the  modified  Neuber  approach) 
along  the  monotonic  stress-strain  curve,  the  surface  layer  will  experience  an  equal  strain  ex- 
cursion along  the  hysteresis  curve  described  in  Appendix  A.  The  local  surface  stress  after 
the  first  loading  is  then  given  by 

<^,nax  = *^0  + 

A A 

where  a (e)  is  the  stress  from  the  hysteresis  curve  at  e.  Subsequent  cyclic  loading 
causes  a strain  excursion  Ae  in  both  the  surface  layer  and  subsurface  region.  The  surface 
mean  stress  is  then  given  by 

a^  = a^+  a (’e)  - Vi  a (Ae)  (B-2) 

The  terms  in  (B-2)  can  be  identified  as  the  initial  stress  state,  the  first  loading  cycle  stress 
excursion  and  one  half  the  cyclic  stress. 

For  the  SC  specimens  which  were  electrochemically  mach  ned,  or  for  an  ideal  bolt  hole  con- 
sisting of  completely  virgin  material,  the  surface  mean  stress  is 

o^=  a (e)  - Vi  o (Ae)  (B-3) 

where  a (e)  is  the  stress  from  the  monotonic  stress-strain  curve  at  e. 

It  should  be  noted  that  o^nax 

®max<®(^  + S>  (B-4) 

f hat  IS.  the  maximum  local  stress  cannot  be  greater  than  that  permitted  by  the  monotonic 
wrrss-sirain  curve.  For  the  materials  considered  in  this  document,  (B4)  can  be  written  ap- 
''  ximately  as 
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(B-5) 


‘'max  ^ ® 

If  (B-5)  is  violated,  e.g.,  for  a particularly  large  first  quarter  cycle  stress  excursion,  the  sur- 
face mean  stress  will  be  given  approximately  by  (B-3).  In  such  a case,  the  effect  of  the  sur- 
face layer  initial  stress  state  is  essentially  “wiped-out”  and  subsequent  behavior  of  the  sur- 
face layer  will  be  identical  to  the  subsurface  behavior. 

Comparison  of  (B-2)  and  (B-3)  reveals  that  the  local  mean  stress  for  a BH  specimen  and  a 
SC  specimen  undergoing  the  same  strain  excursions  differs  by  an  amount 

Aujn  = a (e)  - Op  - a (^ 

Thus,  the  residual  stresses  due  to  machining  effectively  shift  the  mean  stress  values  for  bolt 
holes  relative  to  SC  specimens,  and  it  is  clear  that  this  shift  is  a function  of  the  maximum 
load  applied  to  the  specimen. 


Figure  B-l  Cyclic  Loading  of  a Bolt  Hole  Specimen 
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APPENDIX  C 


MEAN  STRESS  RELAXATION  MODEL 


During  repeated  flights,  a turbine  disk  bolt  hole  may  be  subjected  to  temperatures  and  stresses 
sufficient  to  cause  time  dependent  material  behavior  in  the  hole.  Because  the  bolt  hole  re- 
gion is  constrained  by  surrounding  elastic  material  in  the  disk,  the  hole  region  can  be  viewed 
as  a strain-controlled  region.  In  this  strain-controlled  environment,  the  time  dependent  ma- 
terial response  will  result  in  stress  relaxatioh. 


The  following  procedure  recognizes  the  composite  influences  of  the  various  activities  ex- 
perienced during  repeated  applications  of  a complex  mission.  Any  complex  mission.  Figure 
C-1,  can  be  subdivided  into  separate  activities,  each  of  which  will  contribute  to  the  relaxa- 
tion of  the  bolt  hole  stresses.  Consider  the  stress  history  during  the  i™  activity  of  the  first 
flight  of  the  complex  mission.  The  salient  features  of  the  activity  are  shown  in  Figure  C-2. 
During  the  first  complete  cycle  of  the  activity  (a-b-c-d-e),  the  time  averaged  stress,  is 
given  by 


'lO 


Oltj  + Ut 
ti  +t2 


(C-1) 


where  oj,  02  = the  first  cycle  maximum  and  minimum  stresses  respectively 

*1  ’ ^2  ~ “dwell”  times  at  the  maximum  and  minimum  stresses  respectively 


A cruise  activity  such  as  the  2nd  activity  in  Figure  C-1  can  be  modeled  by  letting  02  = t2  = 0. 


Stress 


Figure  C-1  The  Component  Activities  in  a Typical  Mission 
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Figure  C-2  Simplified  History  of  the  i^^  Mission  Activity 


Equivalently, 


AO: 


(C-2) 


where 


nO  =^ilf2  = 
mi  o 


= arithmetic  mean  stress  during  the  first  cycle 


Ao-aj 


Based  on  the  results  of  the  cyclic  relaxation  tests  reported  in  Appendix  A,  it  will  be  assumed 
that  the  time  averaged  stress  of  the  i event  decays  exponentially  during  the  first  flight. 
That  is 


where 


®il 


Tj  = duration  of  cycling  in  the  i ™ event. 


(C-3) 


Then  the  change  in  both  the  time  average  mean  stress  and  the  arithmetic  mean  stress  after 
time  Tj  is 


Aojj  = 


lO 


(e-“Ti_  1) 


(C-4) 
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Suppose  the  mission  to  be  analyzed  contains  tj,  different  activities.  Then  the  total  change 
in  mean  stress  due  to  those  tj  activities  in  the  first  flight,  <Atr>j,  is  given  by 

V 

<Ao>j  = 2 Aojj 
i=l 

<AOn,>i  = 2 

i=l 

During  the  second  flight  of  the  mission,  the  arithmetic  mean  stress  of  the  first  cycle  of  an 
event  i will  be  shifted  by  an  amount  <Aaj„>j  from  the  corresponding  value  in  the  first 
flight.  Consequently, 


^®i2  = (®io  + <^°m>l)  (e  1) 

V 

<^®m>2  = ^‘'i2 

i=l 

<^‘^m>2  = (®io 


(C-6) 


(C-7) 


In  general,  the  mean  stress  relaxation  during  the  j™  mission  is 


i=l 


<Aa^>j  = 2 (dio  + 2 <Aa^>j)(e-«Ti_,) 

i=  1 k=  1 


(C-8) 


(C-9) 


For  convenience,  let 

A = 2 ffjQ(e-“Ti_  1) 
and  i=l 

^ X. 

B = 2 (e  ■ “ ‘ - 1 ) 
i=l 


--  - - — 
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Consider  now,  the  relaxation  during  the  j and  0-1)  flights 

<Aani>j  = A + B I <Aa^>i  + <Aan,>2  + ' • • + | 

(C-10) 

<^«m>j-l  = A + B { <Aa^>i  + <Aa^>3  + . . . + <Aa^>j.2}  (C-1 1 ) 

Subtracting  (C-1 1)  from  (C-10)  gives  a recursion  formula  for  <Aaj^>j 

<Aa^>j  = (B  + l)<Aaj„>j.l  (C-1 2) 

And  (C-1 2)  can  also  be  written  as 

<Aa^>j  = (B  + l)i-(  <Aan,>,  (C-1 3) 

The  total  change  at  the  end  of  N flights  of  the  mission,  Oj^,  is  then 
N 

0^  = 2 <A0n,>:  (C-1 4) 

j=l 

N 

aR  = <Aa„>l  2 (B  + l)J-(  (C-1 5) 

j=l 

(C-1 5)  is  simply  a geometric  progression  and  can  be  written  as 

(B  + 1)N-1 

aR  = <Aa„>l 1 (C-1 6) 

or 

(^io(e’“'^'-» 

2 (e  - “Ti  _ , ) 

i=l 

This  relaxation  expression  can  be  applied  to  simple  cycle  tests  of  specimens.  For  the  sim- 
ple rapid  cycle  test  shown  in  Figure  C-3a,  (C-1 7)  reduces  to 

(e'"^- •)  (C-(8) 

At  time  T, 


|[  1 +2  (e-"”^!  - 1)  - l! 

i=l  1 


(C-1 7) 
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Thus,  the  mean  stress  of  such  a cycle  would  relax  to  zero,  piven  sufficient  time. 
For  the  dwell  cycle  shown  in  Figure  C-3b,  (C-17)  reduces  to 

At  time  T, 

O 

Thus,  the  mean  stress  relaxes  to  the  limit 


Ao 


and  the  minimum  stress  in  the  cycle,  Omin’  same  amount  to 

o Aa 


Figure  C-3  Cyclic  Test  Characterization 
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However,  it  is  obvious  that  can  never  be  less  than  the  compressive  strength  of  the  ma- 
terial. Thus,  (C-19)  is  valid  so  long  as 


o Ao 

®2  - «m  - “2  ^ °L 


(C-22) 


where  aj^  is  the  compressive  strength  of  the  material. 

Similarly,  (C-17)  is  valid  so  long  as  the  minimum  stress  in  the  flight  does  not  exceed  the 
compressive  limit,  oj^. 

In  summary,  the  relaxation  due  to  N missions  is 

n -r- 

2 aio(e-“ri_i) 

i=l 

s,  1) 

1=1 

for®min>  ^V 


[ 1 + 2:  (e'“'^i  _ 1)  jN  _ 1 
i=l 


(C-23) 
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STATISTICAL  CONSIDERATIONS 
MULTIVARIATE  REGRESSION  ANALYSIS 

The  form  of  the  equation  assumed  to  describe  simple  cycle  LCF  behavior  is 

N = AAe®10^m 

where  N = fatigue  life 

Ae  = strain  range 

= mean  stress 

The  constants  A,  B and  C are  determined  from  a multivariate  regression  analysis.  This  section 
provides  a brief  description  of  such  an  analysis, 

For  each  specimen  test  conducted,  the  fatigue  parameters  Ae  and  are  known  (either 
from  direct  measurement  or  from  analysis)  and  the  resulting  fatigue  life,  N,  (in  this  contract, 
the  number  of  cycles  to  the  appearance  of  a 1/32  inch  crack)  is  observed.  Many  such  tests 
are  conducted,  and  the  assumption  made  is  that  the  typical  behavior  can  be  described  by 
(D-1)  with  the  proper  values  of  A,  B and  C.  These  values  are  most  easily  found  by  consider- 
ing the  logarithmic  form  of  (D- 1 ) 

log  N = log  A + B log  Ae  + Cam  (D-2) 

Equation  (D-2)  is  now  thought  of  as  an  equation  to  predict  the  fatigue  life  of  each  data  point 
from  the  known  loading  conditions,  Ae^  and  aj^^j,  of  that  data  point.  The  difference  between 
the  actual  fatigue  life  Nj  and  this  predicted  life  is  the  error,  Ej,  in  the  prediction. 

Ej  = log  Nj  - log  A - B log  Acj  - Camj 

The  constants  log  A,  B,  and  C are  found  by  minimizing  the  sum  of  the  squares  of  the  errors 
for  the  entire  data  set  (n  data  points),  i.e.,  by  minimizing  the  expression 


2 [IogNj-logA-BlogAej-Ca^.]2  (D-4) 

This  is  done  by  differentiating  (D-4)  partially  with  respect  to  log  A,  B,  and  C,  equating  these 
partial  derivatives  to  zero,  and  solving  three  resulting  equations  for  log  A,  B,  and  C. 

EVALUATION  OF  DATA  CORRELATION 


Fatigue  prediction  systems  quite  frequently  display  a significant  degree  of  scatter  resulting 
from  uncontrolled  or  unrecognized  sources.  Several  of  the  major  sources  of  scatter  can  be 


listed  as  follows: 


1 . Inherent  fatigue  scatter 

2.  Heat-to-heat  material  variations 

3.  Fatigue  model  inadequacies 

4.  Inaccurate  stress  analyses 

5.  Variable  machining  procedures 

6.  Application  to  various  components 

7.  Variations  in  component  usage 


With  so  many  sources  of  variability,  it  is  very  difficult  to  judge  the  contributions  of  a single 
source.  In  this  contract,  these  sources  of  scatter  were  limited  to  as  few  as  possible.  Only 
single  heat  codes  of  material  were  considered:  H/C  CAAZ  in  PWA1216  and  H/C  XNNZ 
in  PWA1057.  Only  disk  bolt  holes  under  Isothermal  loading  conditions  were  considered, 
and  the  variety  of  those  loading  conditions  was  limited  through  extensive  surveys  of  actual 
component  use.  Machining  procedures  for  each  specimen  type  and  the  Ferris  Wheel  disks 
were  controlled  as  nearly  as  possible.  In  addition,  the  most  accurate  state-of-the-art  methods 
of  stress  analysis  were  employed.  Consequently,  the  list  of  major  sources  of  scatter  was 
narrowed  to  the  simple  cycle  and  cumulative  damage  fatigue  models  themselves  (Sections 
V and  VI),  the  surface  layer  and  relaxation  algorithms  (Appendix  B and  C)  used  in  those 
fatigue  models,  and,  of  course,  the  inherent  material  fatigue  scatter. 


It  is  necessary  to  evaluate  the  fidelity  of  the  various  models  and  algorithms  used  to  corre- 
late the  fatigue  data.  A statistical  analysis  of  the  correlation  can  be  used  in  such  an  evalua- 
tion. In  this  contract,  it  is  assumed  that  the  fatigue  data  can  be  characterized  by  a two- 
parameter  WeibuU  distributiuii'^-  The  Weibuil  distribution  is  given  by  the  formula 

F(N)  = l-e'^^/®^  (D-5) 


where  N 

F(N) 


d 


= Random  variable  (life) 

= Cumulative  probability  of  failure  (the  area  under  the  distribution  curve 
from  n = O to  N) 

= Weibuil  slope  (shape  parameter) 

= Characteristic  life  (life  corresponding  to  63.2  percent  cumulative 
failure  rate) 


Figure  D-1  illustrates  two  Weibuil  distributions,  each  having  a characteristic  life  Q = 40,000 
cycles  but  having  different  Weibuil  slopes,  /3,  which  define  the  amount  of  scatter  in  the 
distributions.  For  the  distribution  having  1 .3,  the  BIO  life  (life  at  which  10  percent  of 
the  population  would  have  failed)  is  7500  cycles,  indicating  a fairly  large  degree  of  scatter. 
In  contrast,  the  distribution  with  Weibuil  slope  /)  = 5 has  significantly  less  scatter  since  its 
BIO  life  is  25,000  cycles. 
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Figure  D-]  Sample  Weibull  Distribution 

The  Weibull  probability  plot  then  can  easily  be  used  to  compare  the  overall  scatter  produced 
by  using  a certain  correlation  model  to  the  scatter  that  would  have  resulted  from  a hypothet- 
ical “perfect”  correlation  model.  A perfect  correlation  model  would  predict  the  typical  life 
of  every  test  condition  exactly.  In  such  a model,  the  ratio  of  actualjives  of  specimens  tested 
at  a given  condition  to  the  typical  life  at  that  test  condition  (Nact/N)  would  describe  a 
Weibull  distribution  with  shape  parameter,  /Jp.  All  models  producing  less  than  perfect  corre- 
lation would  have  slopes  smaller  than  when  the  ratio  of  actual  lives  to  lives  obtained 
from  the  correlation  model  (Nact/Nc)  are  plotted.  It  is  desirable  that  P be  as  close  as  possible 
‘o^p- 

Certain  other  features  of  the  Weibull  plot  can  yield  pertinent  information  about  the  correla- 
tion model  being  analyzed.  The  data  points  defining  the  Weibull  line  should  be  continuous: 
sharp  “knees”  or  discontinuous  subsets  of  data  can  signal  poor  correlation  for  certain  data 
subsets.  Likewise,  the  data  should  be  evenly  distributed  on  the  Weibull  line.  If,  for  instance, 
all  dwell  cycle  tests  appear  at  low  failure  rates  while  all  rapid  cycle  tests  appear  at  a high 
failure  rate,  it  would  indicate  that  the  correlation  model  was  too  optimistic  for  the  dwell 
data. 
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DISK  FAILURE  PREDICTION 

The  life  prediction  models  developed  in  this  program  were  based  on  the  prediction  of  typi- 
cal fatigue  behavior.  For  life  prediction  of  an  engine  disk  with  20  bolt  holes,  it  is  necessary 
to  specify  whether  the  prediction  is  to  be  made  for  the  typical  hole  or  for  the  weakest  hole 
in  the  disk.  If  the  prediction  is  to  be  made  for  the  weakest  hole,  then  the  life  predicted  by 
the  typical  model  must  be  adjusted  to  account  for  a failure  rate  corresponding  to  one  in  20. 
This  can  be  done  by  entering  the  Weibull  curve  for  the  correlation  model  at  a median  rank 
corresponding  to  the  first  failure  in  20  holes  (i.e.,  a median  rank  of  3.4  percent'* ) and  find  the 
appropriate  debit  factor,  f,  as  shown  in  Figure  D-2.  The  predicted  life  will  then  be  f x N^,. 


Cumulative 
% failed 


Actual  life 
Predicted  life 

Figure  D-2  Weibull  Plot  Used  to  Determine  Life  Debit  Factor 

It  is  usually  desired  to  predict  the  number  of  cycles  to  failure  (in  the  weakest  bolt  hole)  for 
a large  population  of  disks,  each  containing  tj  bolt  holes.  To  do  this,  an  acceptable  disk  fail- 
ure rate,  Fp,  must  be  specified.  For  example,  it  may  be  desired  to  predict  the  number  of 
cycles  required  to  cause  one  bolt  hole  to  crack  in  10  percent  of  the  disks.  The  life  debit 
factor,  f,  to  be  applied  to  the  mean  life  prediction  is  determined  as  follows; 

The  Weibull  distribution  from  testing  of  single  bolt  holes  is  given  by  equation  D-6. 

F = (D-6) 

where 

F = reliability  rate  in  single  bolt  holes 
f = life  debit  factor  ( = Nact/Nc) 

= shape  factor  and  characteristic  life  of  the  bolt-hole  specimen  Weibull 
plot. 
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APPENDIX  E 


COMPUTER  PROGRAM 

This  computer  program  predicts  component  low  cycle  fatigue  (LCF)  life  exhaustion  as  a 
function  of  user-defined  material  data  base  and  mission  analysis  history  and  includes  the  ef- 
; fects  of  mission  number  and  ordering  of  different  flights.  The  program  embodies  the  results 

of  a research  and  development  effort  described  in  the  main  body  of  this  document. 

I ’ 

The  program  calculates  the  LCF  life  exhaustion  due  to  multiple  flights  of  user  defined  mis- 
sions. Life  exhaustion  due  to  repeated  flights  of  a single  mission  may  be  analyzed  as  well  as 
the  effect  on  life  exhaustion  due  to  a change  in  the  mission  usage.  As  many  as  three  blocks 
of  missions  may  be  analyzed. 

Each  mission  is  described  as  a sequence  of  events  selected  from  a “damage  events  library” 
chosen  to  reflect  various  flight  activities.  Each  damage  event  requires  specific  user  input  (e.g., 
the  maximum  nominal  stress  level)  as  well  as  default  values  of  parameters  which  describe  the 
event  completely. 

Unless  default  parameters  are  over-ridden,  they  will  reflect  the  usage  of  a specific  engine 
disk:  the  TF30-P-1 00  fan  disk  if  a “cold”  analysis  is  specified,  and  the  TF30-P-1 00  third- 
stage  turbine  disk  if  a “hot”  analysis  is  specified.  For  all  events,  nominal  stress  level  input 
must  reflect  appropriate  thermally  and  mechanically  induced  stresses. 

' This  computer  program  is  limited  to  life  predictions  of  isothermal  engine  disk  bolt  holes. 

] Specification  of  a “cold”  analysis  implies  that  stress  relaxation  effects  will  not  be  analyzed. 

, Specification  of  a “hot”  analysis  implies  that  relaxation  effects  will  be  modeled  by  a single 

' , exponential  decay  of  all  stresses.  This  decay  rate  and  a lower  bound  for  the  minimum  allowed 

' I stress,  OL,  must  be  supplied  by  the  user. 

t ! Local  bolt  hole  stresses  and  strains  are  computed  internally  using  a Neuber  calculation. 

Neuber  strain  correction  curves  are  optional  use  ' mput  but  are  recommended  especially  for 
low  strain  hardening  materials.  A surface  layer  algorithm  is  implicit  in  the  calculation  of  bolt- 
' hole  stress.  User  input  must  include  initial  stress  state,  a g.  Component  geometry  is  specified 

through  an  elastic  stress  concentration  factor  and  the  number  of  bolt  holes  in  the  disk  being 
I analyzed.  The  analysis  will  be  performed  for  the  disk  failure  form 

N = AAe®10^m 

I * t 

where  N is  the  typical  simple  cycle  life,  Ae  is  strain  range  and  is  mean  stress.  The  values 
A,  B,  and  C must  be  supplied  along  with  a description  of  the  fatigue  scatter.  The  fatigue 
scatter  is  described  by  the  WeibuU  distribution  parameter,  /?,  i.e.,  the  slope  of  the  Weibull 
distribution  for  specimen  data. 

Double-damage  concepts  are  used  in  the  prediction  of  the  cumulative  damage  effects.  The 
crack  propagation  model  assumes  that  the  material’s  crack  growth  behavior  is  described  by 
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the  standard  Paris  law: 


^ = C (AK) 
dN 


n 


Program  output  is  specific  to  disk  failure  rates  and  mission  ordering  specified  by  the  user. 
A listing  of  the  program  is  presented  on  the  following  pages. 

Warning:  The  data  supplied  in  this  program  is  for  illustrative  purposes  only.  The  program 
not  to  be  used  for  component  design. 


or»r»  ono  non  non 


C DATA  SET  FLMAIN 

COMMON  /lO/  lDl*ICU«MAXtNMAX,MAXCV,lSWAP 
COMMON  /DATA/  CURVM ( 50» 2 » *00 RVHI 50 ,2 ) ,KT , St 21 . 3) t 21 ♦ 3) 

2 fNUMMiS  f SI  GL  CW,  SIGZRO .NEVENT 1 3 » ♦ 6 tNMONO,  NMYST  tNOCC  UR  ( 21  » 3> 

3 ,lH0T,NTPV6t  21,3) 

COMMON  /LOATA/  A,B,C,0,XN,NTlMeSI3  >,ALPHA 
COMMON  /CRKCT/  CRMUNOI  50,2)  ,NCR1,CRHYST|  50,2)  , NCR2 
COMMON  /SRELAX/  ASUR,BSUR tCSUR, ASUB, B SUB , CSUB ,R£L 1 , RE L2 
COMMON  /CEFT/LIFfcSP 
COMMON  /SVLIF  E/XNA,RA,OVA 
COMMON  /ERROR /IEkR 
COMMON  /RSLl/BtTA,FRATE,NBOLT 
REAL  LIFESP 
REAL  KT 

DIMENSION  SURSVEl42,3.4>,ICVCLEI21,3,3),IDAMt3),SLIFE(21 ,3,2 
DIMENSION  I0FF(3) 

READ  INPUT  DATA 

IF(IERR.LI.O)  GO  TO  210 
CALL  INPUT 

INITIALIZE  constants  FOR  STRESS  HISTORY  CALCULATIONS 

1LIN=0 
LIFESP*0.00 
surlst=sigzro 

STRLSTaO.O 
SU8LST=0»0 
JT1M6»0 

NTIMEStNUMMIS  )=0 

COMPUTE  STRESS/S TRAIN  HISTORY 

DO  200  MiSSal tNUMMiS 
00  5 J*l,4 
DO  5 JJ=1,42 
5 SURSVEt JJ,MIS S,J)aO.O 
IFtMISS.EU.l)  GO  TO  20 
STRLSTaSUKSVE  IMJM.MISS-I, 1 ) 

SU6lST»SURSVE IAUM, miss-1, 3) 

SURLSTaSUKSVE  ClUM ,M1SS-1 ,2 ) 

20  CALL  STReSStMlSS,S1RLST,SUBLST,SURSVE,NUM,SURLST,lLlN> 
DETERMINE  STRESS  CYCLES  AND  RELAXATION  FACTORS 
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CALL  CYCLe«SURSVE,ICVCLEtlDAM,IOFF,Ml  SS.NUM) 
C 

C IF  HOf  FLIGHT  RfcLAX  STRESSES 

C 

REL1>0.0 

REL2>C.O 

lFlMlSS.E0«NUMf11S>  GO  TO  35 


30  IFnH0T.E0.2l  CALL  RELAX!  SURSVF, MI  SStNUM  » 

C 

C DO  SIMPLE  LIFE  PREDICTION  AND  PROPAGATION  DAMAGE  DEBITS 

C 

35  CALL  LIFE(SiJRSVE,ICYCLE,IDAM,SLIFE,10FF,MISS) 


C 

C LIFE  EXMAUSTI UN  CALCULATION 

C 

CALL  EXL IFE ! S LIFE, ICVCLE, IDAM, MISS, JTIME ) 

1FIMISS.LT»NUMh1S»  go  TO  200 
JTIME=JTlMEt-l 
C 

C number  of  CYCLES  PREDICTED  NOT  THE  SAME  AS  LIFE  USED  TO 

C RELAX  STRESSES  ITERATION  NEEDED 

C 

DO  50  1=1 ,NUM 

SURSVEU,M15S  ,2)=SURSVE  ( I ,MI  SS  ,2  ) -«E  L 1 
SURSVE(I,MISS  ,3»=SURSVE(1 ,MISS,3)-REL2 
50  CONTINUE 
GO  TO  30 
2C0  CONTINUE 
210  STOP 
END 

SUBROUTINE  INPUT 
COMMON  /ERROR /lERR 

COMMON  /10/I9  I,IDC,MAX,NMAX,MAXCV,ISWAP 

COMMON  /OAMAGE/SDEF!  6,2)  ,R0EF  16,2)  ,TIDEF  < 6,2 ) , T2DEF  16 ,2> 

COMMON  /OAT  A/ CUR  YM(  50,2  ) tCURVHt  50,2>,KT,S  <21,3),R(21,3),  NUMMI  S, 

2 SIGLOW,  SIGZROW^EVENTf3),E,NMOtJO,NMYST,NDCCUR!  21,3) 

3 ,IH0T,NTVPEI  21,3) 

COMMON  /lDATA/  A,B,C,0,XN,NTIMESI3  )♦ alpha 
COMMON  /TIME/TI! 21,3) ,T2( 21,3) 

COMMON  /CRRCT/CRKUNO! 5C,2 )»NCRl,CRHYST<5D ,2), NCR2 
COMMON  /RELI/BcTa, PRATE, NBOLT 
C lERR  -ERROR  FLAG  TO  INDICATE  FATAL  INPUT  ERROR 

C lOl  -LOGICAL  UNIT  FOR  CARD  READER 

C IDO  -logical  UNIT  FOR  LINE  PRINTER 

C MAa  -maximum  number  of  damage  EVENTS 

C NMAX  -MAXIMUM  NUMBER  OF  DATA  POINTS  FOR  S TRE SS /STRA I N CURVE 
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o r-  n o o o o o o n o n o o n r.  n o o r» 


N6V6NT<I 

NTIMES 

MIS 

NTYP6I I.J 
N0CCU«(1, 
SII»J) 

Ri ltJ> 
T1(1,J» 
T2(I,J) 
SOEFIM) 
RDEFiN) 
TIDEFIN) 
T2DEF(N» 
IHOT  = 

I MOT  2 


> -NUMBER  OF  EVENTS  IN  MISSION  1 
-NUMBER  OF  TIMES  MISSION  A IF  FLOWN 
-1  FCR  MISSION  A 

2 FCR  MISSION  B 

> “TYPE  OF  damage  event  FOR  M1SSICN»J 
J)“NUMt!€H  OF  OCCURANCES  OF  EVENT 

-nominal  STRESS 


ITH  EVENT 


“STRESS  RATIO 
“TIME  1 
“TIME  2 

“DEFAULT  VALUE  OF  NOMINAL 
“DEFAULT  VALUE  OF 
“DEFAULT  VALUE  OF 
“DEFAULT  VALUE  OF 
1 COLO  FLIGHT 
HOT  h light 


STRESS  CF  nth  TYPE  EVENT 
stress  RATIO  FOR  NTH  TYPE  EEVENT 
FOR  nth  type  event 
FOR  nth  type  event 


time  1 

TIME  2 


REAL  KT 

DIMENSION  IPRINT(20.3,3),ICHCK(2) tNTITLEi 3,6) 
DIMENSION  IA( B^>,PRINTI3>,t:tLE(80) 

DATA  IBL/lH  / ,1a5K/1H*/ 

DATA  1CHCR/1HM,1HC/ 

DATA  PRINT  /I  HA, IHB, IHC/ 

DATA  NTITlE  /AHTARE,^H  0FF,4H  ,AHTR1M,4H  PaO,AH 

2 AHCRU1,^mSE  .‘♦H  ,4HC0MB,4HAT  ,^H  ,4MT0UC, 

3 AHM  C ,AHGO  ,aHBOMB,AM  RUN, AM  / 

ECHO  INPUT  DATA 

WRITEf:D0,265  » 

1 REAO(1D1,270) lA 

1F1 ECFIIDI) •GTcO.O)  GO  TO  5 
WRITE4IDa,280  ) lA 
WRITE(ISWAP,270) lA 


GO  TO  1 

5  WRITE (100,290  I 
IOI=ISWAP 
REWIND  IDI 
WRITE (IUO,1GCO» 

WRI7E(IDO,1001) 

WRITEf 100,1002) 

WRIT£<1DC,1003) 

WRITE  ( 100  , iu'>A  t 

1000  F0RM4TI IHI,// ,«X,2AHSTATEMENT  OF  LIMITAT  IONS  ,/// ,9X, 

2 A3HTHIS  COMPUTER  PROGRAM  IS  LIMITED  TO  THE  LIF, 

3 52hE  predict  1(XS  OF  ISOTHERMAL  ENGINE  DISK  BOLT  HOLES.  , 

4 19MSPEC1FICAT10N  OF  A ,/,8X, 

5 A3H  COLO  ANALYSIS  IMPLIES  THAT  STRESS  RELAXAT, 

6 33H1DN  EFFECTS  WILL  NOT  BE  ANAL  I ZED . •/,9 X, 

6 1AHSPECIFTCAT10\  , 

7 29h  Of  A »HOT*  ANALYSIS  IMPLIES  • 


o u 


8 ^3HTHA1  relaxation  EFFECTS  WILL  BE  MODELED  BY  ,/,9X, 

9 43HA  SINGLE  EXPONENTIAL  DECAY  OF  ALL  STRESSES.) 

1001  forma T<9X,A^3H MISSION  DESCRIPTION  MUST  BE  SELECTED  FROM  T? 

2 43HHE  EVENTS  LIBRARY.  UNLESS  EVENT  PARAMETERS  ./.RX, 

3 27MARE  OVER-RIDDEN,  THEY  ^ILL  , 

4 45MREFLECT  T HE  USAGE  OF  A SPECIFIC  ENGINE  DISK  :,/,9X, 
t)  3AhTHE  TF30Pi:0  FaN  DISK  IF  A 'COLD*  , 

6 38MANALYS1S  IS  SPECIFIED  AND  THE  TF30100  t/tRA, 

7 43HThIRD  turbine  DISK  IF  A 'HOT*  ANALYSIS  IS  S, 

8 38HPECIFIEO.  rXJMlNAL  STRESS  LEVELS  INPUT  , 

9 34HMUST  REFLECT  APPROPIATE  THERMALLY  ,/,9x, 

I  43MAND  MECHANICALLY  INDUCED  STRESSES.  PROGRAM  , 

1 41H0UTPUT  IS  SPECIFIC  TO  DISK  FAILURE  RATES  , 

2 23HAN0  MISSION  ORDERING  BV,/,9X, 

3 9HTHE  USER.) 

1C02  F0RMAT(//,9X  ,h3HThE  LIFE  EXHAUSTION  MODEL  USED  IN  THIS  PROG, 

1 39HRAM  ASSUMES  THE  FOLLOWING  OESCRI PTION S. ,// , IIX, 

2 5AH1,  typical  SIMPLE  CYCLE  FATIGUE  LIFE, N, IS  DESCRIBED  BY,/, 

3 45X,4CHN=A»< DELTA  STRAIN >**B*10**(C*S1GM A MEAN)  ,/) 

1003  FORMAT!  11X,AIH2.  FATIGUE  SCATTER  IS  DESCRIBED  BY  A /,£I8, 

2 <,5HULL  distribution  IN  TERMS  OF  ITS  SLOPE,  6 ET  A . , 1 1 X , 

3 73H3.  THE  CRACK  PROPAGATION  BEHAVIOR  IS  DESCRIBED  BY  THE  STANDARD. 
^ PARIS  LAi^.,/,  14X,70HSTRESS  RATIO  EFFECTS  ARE  ASSUMED  TO  BE  INCLUO 
BED  IN  Th£  CALCULATION  OF  , / , 1 4X  ,4  IhTH  E ALTERNATING  STRESS  IME.NSIT 
6Y  FACTOR.  ) 

ICC4  FORMAT!/, IIX,  42H4.  MONQTGNIC  AND  HYSTERESIS  STRES 5- STR Al N , 

1 57HCURVES  MUST  RE  INPUT  FOR  THE  SINGLE  TEMPERATURE  AT  WHICH  , 

2 16HTHE  ANALYIS  Wl LL , / , 14X , 13HBE  PERFORME  0.  ,// ,1IX , 

3 57H5.  THE  STRESS  RELAXATION  MODEL  ASSUMES  THAT  ALL  STRESSES  , 

4 52HDECAY  EXPONENTIALLY  WITH  A SINGLE  DECAY  RATE,  A LPHA. , ///, 14X , 

B 51HWARNING  : THE  DATA  SUPPLIED  IN  THIS  PROGRAM  IS  FOR  , 

6 57m1  LLUSTRAT  IVE  PURPOSES  ONLY.  THE  PROGRAM  IS  NOT  TO  BE  USED, 

7 AH  F0R,/,24X ,17hC0MP0NENT  DESIGN.) 

READ!  IDI, 270)  TITLE 
REA0(I0I,350>NUMM1S,IFLAG 

350  FORMAT! 15, 4X, Al) 

IHOT*fl 

IFCIFLAG.EO.ICHCKll) ) 1H0T*2 
1F(IFlAG.E0.1CHCK(2))  IH0T*1 
IrdHOT.GT.O)  GO  TO  4 
IERR=-999 
WRITE(IDO,360  ) 

360  F0RMAT»1X,10<-1H*),37HFATAL  ERROR  EITHER  HOT  OR  COLD  FLIGHT, 

218h  must  be  specified  ) 

4  MISS»0 

1F(NUMMIS.LE.3)  GO  TO  10 
I ERR *-999 
WRITEnDO,300  I 
READ  MISSION  DATA 

10  MISS=MISS+1 

read (101, 100)  F€VENT(MISS>,NTIMES(M1SS) 


o o r» 


IFINEVENTIMIS  SJ-LE.MAXl  GO  TO  20 
IcRR=-‘»9«» 

WR1TEII00,90C  ) NEVENKMISSI  , MISS, MAX 
900  FORMATaX.lOl  1H#),11HFATAL  ERROR, I5,22h  EVENTS  WERE  INPUT  FOR, 15, 
2 8H  MISSION, /,22X,12HA  MAXUMUM  0F,I5,I2H  ARE  ALLOWED) 

20  N=NEVENT<MISS  ) 

DO  22  1J=1,3 
DO  22  lsl,20 

22  IPRINTII,M1SS  ,1J»  = IBL 
00  3C  1 = ]L,N 

READ! 101,110)  NTYPEI I,MISS),DUMMV ,SI I ♦MISS),R  < I, MISS) 

7,T1  II,MlSS),T2fl  ,MISS) 

NN=NTVPEI I,MI SS) 

IFUNMAX-NTVP  Ell, MISS)  )’>NTVPE  1 1 ,MI  SS ) ) 21  ,21,2  5 

21  1 ERR =-999 

WRITEIIDO,910  ) NTYPEI I ,MISS) ,NMAX 

910  F0RMATI1X,10I  1H*),18HFATAL  ERROR  A TYPE,I5,16H  EVENT  WAS  INPUT, 

2 53H  VALID  EVENTS  ARE  GREATER  THAN  0 aNO  LESS  OR  EQUAL  TO, 15) 

25  IFiSf I,MISS}.GT.0.0)  GO  TO  26 
WRITE  1 100,320 )M1$S,I 

26  IFITltl,MlSS) ,UT.0.0)  GO  TO  27 
IFllMUT.EO.2. ANO.NN.EO,4)  GO  TO  29 
TUI  ,MISS)=T1DEFINN,IhOT) 

IPRINTIl ,MlSS,2)slASK 

27  IFIT2lI,MISSH>&T,0.0)  GO  TO  28 
IFUHOT.EQ-2.aND.NN.EO,4>  GO  TO  29 
T2(  1, MISS  )*T2  0£F INN, IHOT) 

IPRINT|1,MISS,3)=IASK 

28  IFIRI I,M1SS).G1.0.0)  GO  TO  129 
IFIIHOT,EO»2»ANO,NN,EO*4)  GO  TO  29 
R 1 1 ,MI SS ) *RDE  F INN, IHOT ) 

1PR1NT|I,M1SS  .Dx^IASK 

GO  TO  129 

29  WRITEfID0,330  ) 
lERR*-999 

129  N0CCURI1,MISS  )=IFIXI  IDUMMY/f  T1  ( I ,MISS  »^T2  1 1 ,MI  SS ) ) ) ♦€  ,5) 
IFINN.EQ.I.OR  .AM.EO.a)  N0CCUR|I,MISS)=1 
IF(NN.EQ.2)  NOCCURI1,MISS)«IFIXIDUMMY40.5  ) 

30  CONTINUE 

IFIMISS.LT.NUMMIS)  GO  TO  10 

CHECK  THAT  THERE  IS  A TAKE  OFF  AND  LANDING  FOR  EACH  MISSION 

00  37  1«1,MISS 
N»NEVENTI I) 

DO  36  J«I,N 

IFINTYPEIJ,!)  .EO.l)  GO  TO  37 
36  CONTINUE 
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WRITE (100,400  )I 

37  CONTINUE 
C 

C ADO  LAAOING  AT  END  OF  EACH  MISSION 

C 

00  38  1=1, MISS 
NEVENTd  )»NEVENTU  »♦! 

J=NEVENT(  I » 

R(J,I )=0.0 
S ( J , I )=0*0 
Tl( J,I)=C.O 
T2( J,I 

NOCCUR( J,I»=0 

38  CONTINUE 
C 

C READ  MATERIAL  DATA 

C 

READ (101,120)  KT, SIGLOW, ALPHA, SIGZRO,E 
AlPHA=ABS(ALPha) 

C 

C READ  RELIABILITT  OaTA 

C 

REA0(IDI,11S>  bETA,FRATE,NBOLr 
REAO(  101,120)  A,b,C,0,XN 
C READ  MONOTONIC  CURVE 

REAO(IOI,1GO)NMOnO,NCR1 
IFJNMCNC.LE.MAXCV)  GO  TO  40 

1 err  *-999 

WR1TE(I0U,920  ) NMOr.O  ,M AXCV 

920  FORMAT)  IX, I0(  IH*), IIHFATAL  ERROR, 15, 27H  POINTS  INPUT  FOR  MONOTONIC 
2 35H  STRESS  STRAIN  CURVE, 12,  8H  ALLOWED) 

40  CONTINUE 

R£A0< 1DI,120)  ((CURVM(I,J) ,J  = 1,2),I=1 ,NMONO) 

IF(NCR1  .NE.O)  REAOdDI  ,120)  ( (CRMONO)  I , J)  , 0=1,  2 ) , 1 =1  ,NCRl  ) 

C READ  hysteresis  CURVE 

REAOdDI  , iC0)NHYST,NCR2 
IF(NhYST.LE.MAXCV)  GO  TO  45 
I ERR =-999 

WRITE(1DD,930  ) NHVST,MAXCV 

930  FCRMATdX  , 1 C TlH* ) , 1 IHFAT AL  ERROR, 15, 26H  POINTS  INPUT  FOR  HYSTERES 
2 22HIS  CURVE, I5,mH  ALLOWED) 

45  CONTINUE 

READ) IDl, 120)  ()CURVH(I,J),J*1,2),1*1,NHVST) 

IF(NCK2«NE.O)  RE ADI IDl ,120) ( ICRHYSTI1,J) , J=1,2),I  = 1 ,NCR2 ) 

C 

C ECHO  MISSION  OATA 

C 

IF) IhDT.EO.I)  WRITE(I00,440)  TITLE 
IF(1H0T.E0.2»  WRITE(1D0,450I  TITLE 
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i WRITE ( ID0«460  )KTtNeOLT 

460  F0RlUTII0X,3HKT>*Ei5.4,/,10X,7HriB0LTSs«I10) 
j NL«NUMMIS-1 

{ IF(NL.EU.O.O)  GU  TO  606 

DO  805  MISS«1 fNL 

MRlTE(ID0»2a3  I FRINT(HlSS>t  NTIHES(MISS> 

WRITE ( 100*430 ) 

! N«NEVENT|MISS )-l 

I 

I DO  800  I«ltN 

KpBNTYPE(I«nISS) 

WRITE f 100*210  1 1* INTITLE (JP»KP),JP«1« 3 >tNOCCUR( ltHlSS)«SI  I «MISS)* 

2 1PRINT(I*HISS*1  )*Rt  I*MISS)*IPRINT(I*MISS*2)*T1(I*MISS)* 

3 IPRINT(I*MISS*3)*T2II*MISS) 

800  CONTINUE 

IF! IHOT.EQ.I)  WRlTEfI0O*38O) 

IF(IH0T.E0.2)  i«ITEa00*390) 

8C5  CONTINUE 

806  N»NEVENTIN1SS  l“l 

WR1TEII00*220  iPRINTiHISS) 
wRITEn00*43C  > 

00  810  Is:l*N 
KPxNTYPEI  I*HI  SSI 

WRITEIIi>0*210)l*INTITLE(JP«KP)*JP>l*3l*N0CCUR(  !*WISS>*Sf  I*MISS)* 

2 IPRIHTII  *MlSS*i>*Rtl*MlSSl  *lPRlNTn  *M1SS  *2)*T1I  I ,M1SS)* 

3 IPPINT(I*MISS*3)*T2n*MISS) 

810  CONTINUE 

IFIlHOT.EO.l)  WRITEfI00*380) 

1F(IH0T,E0.2)  WRITE! 100*390) 

ECHO  MATERIAL  DATA 

WRITE  1100*230  I 

WRITE  1100*240  I SIGLOW*ALPHA*SIGZRO*E*BETA  *FRATE 
E‘E/1000, 

WR1TE(I00*370)A*B«C*D*XN 

WRIYe 1100*250  )IICURVM(I *J)*J«1*2) *1*1 *NMONO) 

IF(NCRl.GT.O)  WRITE! 100*410) ( (CRMDNO I I * J ) * J^l * 2) * I< 1*NCR 1) 

WRITE (IDO* 260  If  I CURVH( I * J ) * J« 1*2) * 1> 1 *NHY ST  I 
IF(NCR2*GT.O)  WH1TE(IOO*420)  ( t CRHYSTI 1*  J I * J«  1 *2)  * I «1  .NCRZ I 
WRITE! 100*470  I 
RETURN 

C FORMATS 

100  F0RMAT(2I5) 
lie  F0RMAT(15*5X*  SEIO.O) 

200  F0RMAT(/*52X*  7HM1SSI0N*1X*A1*1X*11HDESCRI PT10N*/*10X* 

216HMISSZ0N  A 0CCURS*I5*  6H  TIMES*//* ICX. 5 hEVE NT*5X * AMTVPE *4X » 
SIOHOCCURANCES  .2X*7HN0MINAL*8X*6HSTRESS,9X,6HTIME  l*9X,6HTIMe  2*/* 
4«»0X*6MSTRf.SS*  9X*5HRATI0) 

210  F0RMAT(10X,I5*5X*3A4*I3*5X*E10.4*3(5X*A1*  E10.4)) 

220  FORMAT! IHl*/*  52X *7HMI3S ION* IX *Al * IX* 1 IHOE SCRI PTION* //* lOX »5HE VENT  * 
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2 5X»4HTYPE,-»X  ♦10HOCCURANCESt2X*7HN£]MlNALf  fiX,6HSTRESSt9Xf  6HTIM6  1 , 

3 9X,6HTIM£  Z, /,^0X,6MSYRESSf9X,5HRATI0> 

120  FQRMAT(8E10.0  I 

230  rORMAT« iHi,/, iOX,23HMATERIAL  DATA,//! 

2t0  FORMATaOX,19hLOWER  STRESS  B0UN0=,  E15  .A,  1X,3HK  SI,/,  lOX, 

2 6HALPHA=,E1 5-A,iX,5Hl/HRS,/,lOX,15HlNITlAL  STRESS»,£l5 • A, IX ,3HKS 
31 ,/, 10X.2hE=, 2P£15,A,lX,3HRSI,/,10X,5HBeTA»,0PE15*A,/, 

A lOX,  18H01SK  FAILURE  RATE=,F1C*2) 

2'3C  FOKrtAH  IMI,// ,10X,?8HM0NT0N1C  STR  E SS/ STR  A IN  CURVE  ,/ ,10X,  6HSTRESS,9 
2 X,  6HSTRAIN  ,/,llX,3HKSI,12X,5HlN/TN,/,50t5X,2El5.A,/)) 

263  FORMATUHl,// ,i0X,30HHVSTERESlS  STRESS/STRAIN  CURVE  ,/,  lOX  ,6H5TRESS 
2 ,9X,  6HSTRA1N,/,11X,3HKSI,12X,5HIN/IN,/,50(5X,2E15.4,/)> 

265  FORMAT11H1,13HECmO  OF  INPUT,//) 

270  F0RMA7(80Ai) 

280  FORMAT! IX ,80A I) 

320  FORMAT!  IX, 101  ln*l  ,36HFATAL  ERROR  NOMINAL  STRESS  LEVEL  NOT, 

2 32H  INPUT  FOR  MISSION, 15, IX, 5HEVENT, 15) 

330  FORMAT! IX ,10! IH») ,^5HFATAL  ERROR  DAMAGE  EVENT  TYPE  A HAS  NO  OEFAUL 
2 21hT  values  FCA<  HOT  RUN.) 

29*!  FORMAT!  IX, 12H  END  OF  ECHO) 

30:  F0RmaT!1X,1C.!  Ih*),38HFATAL  error  MORE  THAN  3 MISSIONS  INPUT) 

380  format  !AOX,37H»OEFAULT  VALUE  FOR  TF30-P100  FAN  DISK) 

390  FORMAT! AO X,A1h»0E FAULT  VALUE  FOR  TF30-P100  TURBINE  DISK) 

^00  FORMAT!  IX,  10!  Ih*),35hFATAL  ERROR  NO  TAKE  OFF  FOR  MISSION, 15) 

910  FGRMAT!lHi,/, 10X,36MCORRECTION  CURVE  FOR  MONOTONIC  CURVE, 

,1  /,10X,6mSTRA1N,9X,5HCFACT, 

1 /,10X,5H1N/IN, 

2 /, 501  5X,  2E15  •*♦,/) ) 

920  FORMAT!  IHI,/,  10X,36HCORRECTION  CURVE  FOR  HYSTEP.SIS  CURVE, 

1 /,ICX,6HSTRA  l6i,9X,5HCFACT, 

2 /,10X,5H1N/IN, 

3 /,50I5X,2£15 .9,/) ) 

115  F0RMAT!2E10.0,I5) 

930  F0RMAT190X,3HKSI ,27X,3HHRS, 12X,3MHRS/ ) 

AAC  FORMAT! IHI,// ,52X,21HC0LD  MISSION  ANA LYS 1 S. /, 2 5X , SOAl , // ) 

95C  FORMATdHl,// ,62X,20HH0T  MISSION  ANALYSIS  ,/,25X,8CA  1,//) 

370  FORMATllHl,// ,10X,23HL1FE  EQUATION  CONSTANTS,// 

2 ,10X  ,2HA  = ,E1  5,A,/,10X,2MB«,E15.A,/,10X,2K»,  E15.A, 

3 ////,10X,22HCRACK  GROWTH  CONSTANTS,//, 10X,2HD*, 

9 E15.A,/,10X,  2»-*»*,E15.A) 

970  FURMATi 1H1,17X,27HRESULTS  OF  MISSION  ANAL YSI S , //, 1 0 X, 7HM I SS lON.BX 
2, 10hCUMULAT1VE,6X ,9hNUM8ER  OF ,/ , 2 5X , 1 CHPERCENT  OF  , 

3  6X,7nFLlGrtTS,/,25X,9HLlF6  USED,/) 

ENO 


L 
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SUBROUTINE  ST  RESSt  MISS  t STRLST  ,SUBL ST,  SU» S Vc  ,NUKLST  , SU«»LS T ♦ I L I N ) 
COMMON  /10/I0I,lU0 

COMMON  /UATA/  CURVM( 50,2),CURVH(50,2 ) ,KT, SS<21 ,3) ,R R( 2 1, 3 ) , NUMMI S 
2,SIGLCV,SIG?R  0,NevENTri ) ,E ,NMONO,NHVST  ,NOCCUR  ( 21 , 3> 

3tlHOT,NTYJ>E(2  1,3  ) 

COMMON  /CRRCT/  CKMQNO( 50,2) ,NCR1,CRHYST< 50,2» ,NCP2 
REAL  RT 

DIMENSION  SURSVE  (42,3,4) 

THIS  ROUTINE  CCMPUTES  THE  STRAIN  HISTORY  OP  A MISSION 
INITIALIZE  CONSTANTS 

ILIN=0  STRESS  LEVEL  ON  LINEAR  PORTION  OF  HYSTRESIS  CURVE 
ILIN»1  ON  NON-LINEAR  PORTION 

F^UMsl 

RWIPE=0 

1=1 

S3=0.C 

S^*0.'T 

ISIGN=-) 

ITIME=0 

SURSVEi VUM,MI SS,1 ?=STRLST 
SURSVEIN«JM,rtISS,3)=SU8LST 
SIGSUR»SURLST 
SURSVE(NUH,MI SS,2)=S1GSUR 
1 S«SS(1,M1SS) 

R=RR(l,MISSl 

lCrtNGE*0 

1F(S.gT,S4.AND,1SIGN,LT.0>  ICHNGE*! 

IF(S.LT,S4.AND.IS1GN.ST.0)  ICHNGE  = 1 
IF(ILlN«Eb«0.  CK.ICHNGE.EQ.O)  GO  TO  200 
NUMSCK=NUM 

I FI  I TIME .EO. I .AN0*NUM,E0.1)  NUM*NUMLST 
STRLST*SURSVE IN0M,M1SS,I ) 

SU«lST»SURSVE |\UM,MISS,2) 

SU8lST»SURSV6  <NUM,MISS,3) 

NUM»NtWttCK 

S3»S4 

IL1N»0 

200  CONTINUE 

IF  PART  MAS  NO  PREVIOUS  STRESS  HISTORY  SUBSURFACE  FCLLCMS 
MONOTONIC  STRESS  STRAIN  CURVE 
IF(  I.EO.I.ANO.MISS.EO.I.ANO.ITIME.EQ.OGO  TO  5 
GO  TO  10 

FOLLOW  MONOTONIC  CURVE 
5 DV0£«IRT*S)*( aT*S)/E 

CALL  INTERlCUPVH,NMONOfDE,OVOE,CRMaNO,NCR  I) 

1L1N«1 

SIGSUB«OVOE/OE 

STRAIN«0E 


NUM*NUMLST 


I 


ISIGN=1 
IMONO=C 
lSIGr>j  = l 
GC  Kj  ?0 

C follow  HYSTERtSIS  CURVE  FOR  SUBSURFACE  STRESS  LEVEL 

IC  OVUE=«M*(S3-5.))  »*2/£ 

I SIGN=1 

IF(S3,GT.S)  IS1GN=-1 

CALL  IMTER<C'JkvH,NHYST,DE,DVDE,CRHVST,NC«  ?) 

1LIN=0 

It-<DE.GT.CURVH(i,l))  1LIN=1 
USlGMA=DVDE/Ot 
SlGSue=SUBLST  ♦ISIGN»DS IGM A 
STRAIN  = STRi.ST  ■*^ISiGk*OE 
1FISIGSUB.lt. S1&L0«)SIGSUB=S1GL0W 
C CHECK  THAT  MONCTONIC  CURVE  IS  NOT  VIOLATED 

lMONO  = 'i 

CALL  vLATE«CUKVM,STRAIN,S1GSUB,1VIOL  *NMONC) 
IFHVIOL.lT.O  > go  TO  5 

C HAS  surface  layer  been  previously  BEEN  DESTROYED 
20  IFtKWiPE.Nc.O  ) GU  TO  45 
IFUMCNO.GT.O  It  GO  TO  30 

C SUBSURFACE  FCLLOwED  MOnCTQNIC  STR ESS/STR A IN  THEREFORE 

C IS  SIG  2R0  PLUS  SIG  HYST 

CALL  INTERYI CURVM.NHYST, STRAIN, SIGSUR ) 

SIGSUK=SIGSUK  ♦SIGZRO 
GO  TO  40 

30  S1GSUR=SURLST  ♦ISIGNYDSIGMA 

I F<  SIGSUR. LT. SIGLOw)  S I6SUR  = S IGLOW 
C CHECK  that  MONOTuNIC  CURVE  IS  NOT  VIOLATED 

40  call  VLATc  (CURVM, strain, SIGSUR, IVIOL  ,NM0ND) 

IF( IVIOL.GE.O » GU  TO  50 
K^IPE  = 1 

45  S1GSUR=S1GSUB 
50  IF(R.LE.O)  GO  TO  60 

t)VDE  = (KT*SA'«  1 .0-R))**2/E 

CALL  1NTER(CUKVH,NHYST,DE,0VDE,CRHYST,NCR  2) 
iLlNxO 

IF(0E,GT.CURVhI1,1))  ILIN=1 

DV=OVnE/DE 

NUHxMUK^l 

SURSVE{NUM,MI SS,1)=STRAIN 
SUKSVE(NUM,H1 SS,2>=SIGSUR 
SURSVE<NUM,MI  SS,3)=S1GSUB 
SIGSUR*SIGSUR-OV 

IF( SIGSUR.lt. SIGLOW)  S IGSUR=S IGLOW 
SIGSUBxSIGSUB-OV 

1FISIGSUB.lt. SIGLOW)  SIGSUB= SIGLOW 
STRAIN=STRAIN-Dt 
60  NUM*NUM+l 

SURSVEINUM,MI SS.l I^STRAIN 
SURSVE<'AJM,MI  SS,2)=S1GSUR 
SURSVE*N0H,?AI  S£,3)=SIGSUB 
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SURSV6(NUM,M£  SS,<*»=FLOAHWCCUR(  I ,MlSSn 
IFJITlMb.EU.3  I S4=:i 
1=I  + 1 

1F( iTIME.Eg,! » GO  TO  90 
IF<  I.lE.NEVEMTIMISS)  > GO  TO  1 
1T1ME*1 
I*l 

r4UMLST»NUM 
1VMAX=1 
70  I=I^1 

IFd.GT.NUM)  GO  TO  80 

IF(SURSVE(I,M ISS,2).GT.SURSVEI1VMAX,MISS, 2) ) IVMAX*1 
GO  TO  70 
80  1 = 1 
NUM=1 

90  1 FfNUM.LT.lVM AX)  GO  70  I 
RETURN 
END 

SUBROUTINE  IN  TER  ( ARRA  V ,N,  XB  , X VB  ,CRCT  » ICRC  T) 

COMMON  /lO/  1 Dl,IDO 

DIMENSION  ARRATiSOtZ) t xy( 50) »CRCTI50»2> 

THIS  ROUTINE  FINOS  ThE  VALUE  OF  XB  GIVEN  XVB  SUCH  THAT  THE  POINT 
IS  ON  THE  CURVE  DEFINED  BY  X*Y  VS.  X 

COMPUTE  XY  ARRAY 

DO  10  I=lfN 

10  Xvn  )»AVRAY)  I ,l)*APPAYt  1,2) 

DETERMINE  IF  INTERPOLATION  OP  EXTR APOLATI  CX-t  IS  NEEOFO 

IF(XY»N)-XYS)  30, 2C  ,40 
20  XB*ARRAY< N,1 ) 

GO  TO  909 
C 

C EXTRAPOLATE  SOLUTION 

C 

30  NN»N 

GO  TO  70 
C 

C INTERPOLATE  SOLUTION 

C 

-0  DO  50  I«1 ,N 
NMIGH«1 

IFfXYB-XYlD)  60,60,50 
50  continue 
GO  TO  30 
60  NN«NH1GH 

70  IFINN.EO*!)  G C TO  75 

XM«UPRAV|NN,  2)-ARRAVINN-l,2)  »/l  ARRAYINN,  1) -ARRAY  I NN-1  ,I  ) ) 
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c 


c 

c 

c 
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c 

c 
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c 

c 
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B=A«(<AYINN-lf  ARRAY«NN“1,1)  ) 

GO  TO  80 

75  XM=AR«AYfNM»2  )/ARRAYfNN,l » 

B=0.0 

80  XB=( SQRT(B*B  + 4oO*XM»XYB  >-B ) / f 2.0*XM) 


USE  CORRECTION  CURVE  TO  DETERMINE  ACTUAL  STRAIN  FROM  NUE0ER 

strain  prediction 


999  CFACT=1.0 

IFUCRCT.GT.O  J CALL  1NTERY< CRCT, ICRCT  ,XB  t CFACT  ) 
X8=X3»CFACT 

CALL  INTtRYI ARRAYtN.XBtY) 

XY6=XB*Y 

RETURN 

END 

SUBROUTINE  IN  TERY (CURVE ,N ,X tY ) 

DIMENSION  CUR  V£| 50,2) 

COMMON  /lO/  1 D1,1D0 

This  routine  finds  y given  x using  curve  (1,i)  = xd) 
AND  curve  (1,2  > = Yd  ) 

determine  if  INTERPOLATION  OR  EXTRAPOLATION  IS  NEEDED 
NTOP=N 

IF)CURVE( N,l) -X)  40,  10,  20 
1C  YsCURVECNTUP,2) 

GU  TO  999 

INTERPOLATE 


20  DO  25  1=1, N 

NTOP=I  j 

1F(X-CURVE(I,  1 »l  30,  10.  25 

25  CONTINUE  ! 

GO  TO  40  j 

30  IHNTOP.EO.l)  GO  TO  35 

Y»(  1 CURVE  (NTO  P,2  l-CURVE  )NT0P-1,2  > )/( CURVE  (NlOP  , 1 > -C  UR  Vf  ( NTOP-1 , 1 > ) i 

2)*( X“CURVE(NTCP-1, 1 )) ♦CURVE (NTOP-1  ,2  ) { 

GO  TO  999  ^ 

! 

ASSUME  CURVE  STARTS  AT  0.0, 0,0 


35  Y:*(CURVE(  1,2)/CURVE(1,1))*X 
GO  TO  999 

EXTRAPOLATE  TC  SOLUTION 


*.0 

999 


NTOP=N 
GO  TO  30 
CONTINUE 
RETURN 
END 


j 
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BLOCK  04TA 

BLOCK  DATA  ROUTINE  TO  SET  DEFAULT  VALUES  OF  I/O  DEVICES  ANO  DAMAGE 
EVENT  default  values 
COMMON  /ERROR /lEftR 

COMMON  /lO/  I D1«100«MAX«NHAX*MAXCV«ISWAP 

COMMON  /damage/  SDEFI  6,2),RDEF|  6*2)tTlDEF(  6,2»,T20EF(  6,2) 

Data  101,100  /B«6/ 

DATA  S0EF/i2*0.0/ 

DATA  TlOEF/,3  33,  "i.  Cl  66  , .333 ,3*0.0 166  , .033  3,  ,00833  ,.  5,  .0,  .08333,.  05 
2 / 

0 AT  A RDEF  /O.  . '.5,0*, *83,  *29,  *36, 0*,  .56,0*,  0*,*16,  .45  / 

D AT  A T20E  F /O.  , .0166 , 0 . , .0 16  6,. 0166,  .0 166 , 0.  , . 0 166  ,C.,0..  • 08  333,.  05 
2 / 

DATA  MAX,IERR  ,NMAX,MAXCV,ISWAP 
2 /20, a, 6, 50,1/ 

END 

SUBROUTINE  VL  ATE  I CURVE , STRA IN ,S1GSUB , IVIOL,NCURV 1 

THIS  ROUTINE  DETERMINES  IF  GIVEN  STRAIN  £ SIGSUR 
IF  CURVE  IS  violated 
common  /lO/  1 Ol,lDO,MAX,NMAX,MAXCV 
DIMENSION  CURVE!  50,2) 

XPRT«FLOaT  tNC  URV ) 

IVI0L=0 

CALL  INTERVIC  URVE ,NCURV , STRAIN, STRESS > 

IFISIGSUB.GT. STRESS)  IV10L«-999 
XPRT»FlOaT!  IV  lt!L) 

RcTUKN 

END 

SUBKOUT IN  £ CV  CLE 1 SURSVE ,1 CYCLE , 10  AM, 1 OFF , hI SS ,NUM ) 

DIMENSION  SUKSVE44?,3,4),ICYCLE121,3,3),I0AM(3 ) 

DIMENSION  10FFI3) 

This  subroutine  determines  stress  cycles  and  number  of  times 
That  each  is  cycled 
♦♦BASED  ON  subsurface  STRESS  CYCLES 

ALSO  CALCULATED  ARE  CONSTANTS  FOR  STRESS  RELAXATION 

ICYCLEI I ,M1SS ,J) 

J*l  POINTER  TU  SURSVE  FOR  BEC.  STRESS 
J»2  POINTER  TO  SURSVE  FOR  END  STRESS 
J=3  NUMBER  OF  CYCLES 

lOAM(MISS)  COUNTS  OF  STRESS  CYCLES  IN  A MISSION 

COMMON  /LDaTA/  A,e,C,0,XN,NTlMESI  3),ALPHA 
COMMON  /IQ/  in,  I DO 

COMMON  /DATA/  CURVM 1 50 ,2 ) ,CURVHI 50,2 ) ,KT,  SI 21 , 3 ) , R ( 21 , 3) ,NUMM IS 

2 ,SIGLCW,  S1G2  Rt;,NEVENT  I 3)  ,E,NM0N0,NHYS7,N0CCUR  12)  ,3  ) 

3 ,IHOT,NTYP66 21,3) 

COMMON  /SRELAX/  ASUR , BSUR ,CSUR ,ASUB, B SUB  , CSUB , REL 1 , RE L 2 
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COMMON  /TlM£/TrH21,3>  tTT2<21,3) 

REAL  KT 
ASUR^O.O 
BSUR=0»0 
CSUR  = i.O 
A5UB=0,0 
eSU8=0.0 
CSUB=l.O 
10  isrR=o 
iEveNr=o 
idamimiss)=o 
20  ISTR=1STR+1 

IFIIEVENT  + I.GE.NEVENUMISS))  GO  TO  999 
C 

C IS  THIS  THE  END  OF  A DAMAGE  EVENT 

C 

lF{SURSVEnSTR»MlSS,4)  .EQ.O.O)  GO  TO  20 

C 

C INCREMENT  DAMAGE  C^CLE  COUNTER 

C 

I£VENT  = 1EVENT 
T isTTineVENT  ,MISS> 

T2=Tr2< lEVENT  ,MISS) 

TRAT=tT2-Tl»/ <ll^T2> 

TI  = IT1+T2  »»SUttSVE(ISTR,MlSS,4) 

ESTUFF  = < 1.0/E  XP(  ALPHA *Tin -I  ,0 
C 

C determine  EVENT  TYPE  AND  CALCULATE  NUMBER  CF  CYCLES 
C 

N-NTVPEHEVEN  T.MISS) 

IFIN.NE.A-.AnO  .N.N£.5;aND.N«NE.6)  GO  TO  30 
r 

C GET  NUMBER  OF  CYCLES  DIRECTLY  FROM  SURSVE 

C 

25  SIGM=f SURSVE4  1STR-! ,MI SSt 2 > ♦SURSVE » I SIR. M 1 SS. 2 M /2 . 
DSIG=ABS<  SURS  VEI  I STR-1  f Ml  SS  t 2 )-SUR SVE  US TR.Ml  SS,  2 ) } 
SIGI=SIGM-|OS I&/2.0»»TRAT 
ASUR  = ASuR4SIG I ye ST  off 
BSUR*BSUft^E51  OFF 
CSUR*CSUR  + ESTUFF 

S1GM=« SURSVE < lSTK-1, MISS. 3) ♦SURSVE {ISTR.M ISS. 3 ) )/2, C 

DSIG*ABS( SURS  VE( lSTR-1 ,MISS.3 )-SURSVE ( 1ST R ,MI SS,3 ) ) 

SIG1=S1GM-I0S IG/2.0)*TRAT 

ASUB«ASUB^S1G  li*tSTuFF 

BSUB-BSUa^EST  OFF 

CSU8*CSUB^EST0Ff 

lF«N.fc0.2)  GO  TO  35 

IDAMJMISS  )*10  AMIMISSMI 

I>IDAMIMI SS; 

ICYCL Ell. MISS  .1>=ISTR-1 
ICYClEII.MISS  .?»*1STR 

I CYCLE  1 1.  MI  SS. 3)  =I  F IX  I SURSVE  1 1 STR, MISS, 4)  ) 

GO  TO  20 
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30 


IFIN.Nb.Z)  GO  TO  4C 
GO  TO  25 
35  CONTINUE 

IOA(«(hISS)sIO  AH(MlSSt-»l 
I*I0AMMIS5) 

ICYCLEI  l.MISS  •DsISTR-l 
ICYCuEUtHlSS  .2I=ISTR 
ICYClEU.MISS  ,3>  = l 
GO  TO  20 

IF<N«NE*1  ) GO  TO  60 

FIND  VHAX  and  VMIN  OF  MISSION 

IMAX=2 

IMIN=2 

ICT»2 

45  ICT=ICT-H 

IF<  ICT.GT.NOM  ) GO  TO  50 

IF<SURSVE  fICT ,MISS,2).LT,SURSVE<IMIN,MISS,2»)  IMIN*  ICT 
IFiSuRSve lICT  ,rtISS»2>.GT.SURSVE(IMAX,MlSS  t2»>  IMAX=ICT 
GO  TO  45 

53  I0AM(MlSS)sI0AMfMlSS>4l 
I=>I0AM(MISSI 
I0FF<MISS)»I 
ICYCLEdfMISS  ,1)=IMIN 
ICYCtE(l»MISS  ,2>sinAX 
ICYCLEU  tMISS  ,3}-l 

S1GM=(SURSVE(  lMAX,MlSSt2>-»SURSVE(  iMlNyMI  SSt2)  )/2.0 

OSIG®ABSI  SURS  Vtl IMAXf MISS 1 2 ) -SURSVEf I MIN ♦ Ml SS  t2) > 

SIGI*S1GM-(0S1G/2*0)*TRAT 

ASUR  = ASUR*SI&l»£STUFF 

BSUR*BSUR^ESTUFF 

CSUR=CSUR4ESToFF 

SIGMs(SURSVEI 1MAX,m1SS*3>4SURSVE( IMIN  »MlSSt3) )/2.0 

DSIGsABSI SURS  YE( IMAX «MISS *3 ) 'SURSVE ( I MIN » Ml SS • 3) ) 

SIGIxSIGM-|0S1G/2.0)4TRAT 

ASUB«ASoB4S16 IxESTUFF 

8SUB=BSUB4ESTUFF 

CSUBxCSUB+ESTUFF 

GO  TO  2'» 

6D  IFfN,NE.3)  GO  TO  20 

IFUEvENT  .EQ.D  GO  TO  20 
NL*NTYPE I lEVE NT-1, MISS) 

lF(NL,NE.l.AND*r:L.NE«4.AN0*NL»NE.5.AND.NL  .NE.6)  GO  TO  20 
NN=NTYPE»  IEVENT*1,M1SS) 

lFfNN.NE»l.ANU.NN.NE«4»AND.NN.NE.5.AND.NN  .NE.6  ) GO  TO  20 

1 0AM (MISS )xIO AM(MlSS)4l 

IxTOAMIMlSSI 

ICYCLEI l.MISS  ,1»*ISTR 

ICYCLE< l.MISS .2(=ISTR41 

ICYClE<I,MISS,2»=1 
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SIGM=Sl»9iiVE<  I MR  .MISS, 2) 

SIG1=SIGM 

ASUR  = ASL'K +516  1*6  STUFF 
8 SUR  = B SURGES!  UFF 
CSUR=CSUR+ESTUFF 
S1GM=SURSVE(I  SIR, MISS, 3) 

SIGlsSIGM 

ASUB=ASUB+S1G I*ESTUFF 
0SU8*8SU6+EST  OFF 
CSUB=CSUB+6STUFF 
GO  TO  20 
999  continue 
RETURN 
END 

subroutine  LIFE! SURSVE , ICVCLE .IDAM.SL IFE , lOFF ,M1SS» 

this  SUBROUTiNt  DOES  SIMPLE  CYCLE  LIFE  CALCULATION 
USING  SURFACR  stress  HISORY 


SLIF6  n.MlSS,  J I - LIFE  PREDICTION  FOR  ITH  STRESS  CYCLE  OF  MISSION 
SLIFEU.MISS,  2)  - PROPAGATION  DAMAGE  DEBIT 

IOFF{MISS)  - LOCATION  OF  TAKE-OFF  CYCLE  POINTERS  IN  CYCLE  ARRAY 

COMMON  /lO/  101,100 
COMMON  /SVLIFE/  XNA.RA.DVA 

COMMON  /DATA/  CUKVM(50,2),CURVH| 50,2 ) ,KT, SI2I ,3),RI 21 ,3) .NUMMIS 

2 ,SIGlOW,SIGZRO,NEV6NT<3) ,E,NM0N0,NHYST,N0CCUR  »2I,3  ) 

3 ,IM0T,NTYPE( 21, 3» 

COMMON  /LDATA/  A ,B  ,C  ,0  , XN  ,NT1  MES  < 3 > 

REAL  KT 

dimension  I0FF(3  ),SURSVE  (42, 3,4  >,  I CYCLE!  21,3,3  > 

DIMENSION  I0AH(3»,SLIFEI2I,3,2) 

A,B,C  ARE  USER  INPUT  CONSTANTS  USED  IN  LIFE  CALCULATIONS 

LOG(N>=lOG( A »^B*L0G(0ELT  STRA1N>*C*STRESS  MEAN  WHERE  N*L1FE 


DETERMINE  HOW  MANY  STRESS  CYCLES  IN  MISSION 
N*1DAM!M1 SSI 

LOOP  THROUGH  STRESS  CYCLES  CALCULATING  LIFE 

DO  100  1*1, N 

II  = 1CYCLE(I,MISS,1) 

I2=ICYClE(1,M1SS  ,2) 

DSTR*ABS(SURSVE(I1,MISS,I  )-SURSVE  ( 12  , MISS  , 1 H 
SIGMx(bUKSVE(  II, MISS, 2)*SURSVE(I2, MISS, 2>  )/2. 
ARGxALOGl  01 A ) ♦b*  ALOGlOf  DSTR ) •»C«SI  GM 
SLIFEd.MlSS,  1»*10.**ARG 
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It=(SLlFE<I,HI  SS«I).GE.1.0E6)  SLIF E ( I « MISS  • 1 )>=  1 •0E6 
100  CONTINUE 


CALCULATE  PROPAGATION  DAMAGE  DEBIT 


IF(M1SS»GT.1)  GO  TO  125 
ISUB  = IOFFJrtIS5> 

XNAsSLIFE  nSUBtMlSStl) 

Jl«ICYCLE(ISUBtMlSS«l) 

j2sICyCLE(ISUE>*MlSS»2) 

RA>AMlNlf SURS VE|Jl*MISSt3  >. SURSVB I j2  ,MISS  t3)l/ 

2 AMAXKSURSVE  ( J1»M1SS«3)«SURSVEIJ2«M1SS«3  )) 

DVA«SURSVEIJ1  *HISS«3)->SURSVE(  J2*MISSf  3) 

DVAsABSIDVA) 

125  CONTINUE 

DO  150  lsl»N 
I IsICvCLEll.MiSS,!) 

I2*ICVCLE<ItM lSSt2) 

XNB«SL1FE<I«MISS»11 

RBsAMlNKSURSVEf  lltMlS$f3)«SURSVE(I2«MlSS  »3M/ 

2 AMAXIISURSVE lll*MlSS«3>«SURSVEfl2tMlSSt3>) 

DVB«SURSVEiIl tMISS«3)-SURSVE(I2«MISS*3) 

0VB*ABS(0VB) 

SLlFEd.HlSS*  2>s|  t HO-RB)*OVA  )/l  (D-RA  )*DV  bl  )**XN)*<  XNA/XNB) 
150  CONTINUE 
«»9St  CONTINUE 
RETURN 
END 


SUBROUTINE  EX  LIFE ( SLIFE tICVCLE*IOAMtM ISS « JTIME > 

COMMON  /LEFT/  LIFESP 

COMMON  /DATA/  CURVMI  50*2  ) tCURyHI  50«2  1 «XKT  ,S  1 2 1 «3 ) .R  12 1 *3  ) tNUMMlS  • 

2 SIGLON«SIGZRO,r4EVENT(3),EfNMONO,NHrST,NOCCUR  ( 21«B>  t 

3 IM0T,NTYPEI2 1.3) 

COMMON  /lO/IDl.IDO 

COMMON  /LDaTA/  A.B.C.O.XN.NTIMESIB) 

COMMON  /RELI/  BETA.FRATE .NBOLT 
REAL  LIFESP 

DIMENSION  SLIFE(21.3.2),ICVCLE(21.3.3)«IDaM(3)  ,PRINT(3) 

DATA  PRINT  /I  HA.  IHB.  IHC/ 

NDAM«I0AM (MISS) 

XB0LT«FL OAT (NBOLT ) 

R0»1.0-FRATE 

XNK«( 1.0/BETA  )*(ALOG( ALOG( l.O/RO) )-ALOG(  XBOLT ) 

2 -ALOG( AL0G(2 .0) )) 

XNRsEAP(XNR) 

IF(MISS.EO.NUMMIS)  GO  TO  200 

THIS  IS  NOT  THE  LAST  MISSION  CALCULATE  REMAINING  LIFE 


lOS 


non 


! 


I XLIFE=0.0 

DO  20  1=1»N0A« 

XLlFE  = XLlFE+<  FL0aT<1CYCLE(I  ,MISSt3))  / (SLI  FEd^MlSS,  1»>* 

I 2SLIFEI!,MISS,  2)  » J^NTlMESI  "^ISS  )/XNR 

1 2C  CONTINUE 

IFIXLlFt+LIFE SP-1.0)30,40»50 
, 30  L1F£SP=XLIFE+lIFESP 

' XP«  r*ILlFES»») 

I WRITE (iDOtlOO  IP«INT(MlSS»tXPRT,NTlMES (MISS) 

( RETURN 

' 40  WRITE(IDO,100  )PRINT(MISS),XPRT,NT1MES(MISS) 

WRITEdDOtllO  ) 

110  F0RMATdX,23MLlFE  HAS  BEEN  EXHAUSTED.  ) 

GO  TO  999 

LIFE  USED  UP  BEFORE  MISSION  OVER  CALCULATE  APPROXIMATE  LIFE 

50  XLIFE=IXLIFE)/NTIMES(MISS) 

N = ( 1.0-LlFESP  )/XLlFE 
WRITE (IDO, 120  IPRINTCMISS) ,N 

120  F0RMAT(1X,49HLIF6  EXHAUSTED  BEFORE  SPECIFIED  NUMBER  OF  MISSION, IX, 
2 Al.lX.lHMWERC  COMPLET£Dt/tlXfl3HAPPROXlMATELY,I5tl  3m  MISSIONS  COU 
232HLD  BE  FLOWN,  REANALYSIS  ADVISED.) 

GO  TO  999 

r 
v 

C PREDICT  HOW  MANY  MISSIONS  CAN  BE  FLOWN 

C 

200  XLIF£=0.0 

DO  210  1 = 1, NO  AM 

XLIFE=XL  !»-£♦(  FLCATdCYCLE(I,MISS,3))/(SLI  FE  d , MIS  S,  ? ) ^SL  I FE  ( 2 ,M1  S S 
2,2)) )/XNR 
210  CONTINUE 

N*(  1.0-LIFESP  )/XLIFE 
XPR7=100. 

IFdHOT.EO.l)  GO  TO  220 

IF(NT1MES(M1S  S).EO.O.AND.  JTIME.EO.O)  f-0  TO  950 
NT=NTImES(MIS  S)-N 
IFInT.GT.SO)  GU  to  950 
220  N»(N^5)/10 
N=N*IG 

WRITE  dOn,  100  )PRINT(  MISS)  tXPR  T,n 
GO  TO  999 

P5C  NTIMES»M1SS)=N 
RETURN 

IOC  F0RMATd3X,Al  ,11X,F8.2,7X,I10) 

099  STOP 
END 

♦ 

! 
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i SUBROUTINE  RE  LAX t SURSVE ,Ml SS t NUM ) 

COMMON  /SR£LAX/ASUK*BSUR*CSURtASUB«BSUByCSUB,RELl»RI12 
COMMON  /LOATA/  A *6  ,0 f D , XN  ,NT1  ME S ( 3 ) t A LPHA 

COMMON  /DATA/  CURVM(50,2>.CURVH<  5C,2  » »X<T,S<21  .3)  tR  (21.3  ) tNUMMIS, 
! 2 S1GLOW,SIGZRO.NEVENT(3),E.NMONO,NHYST.NOCCUR(21,3)  ♦ 

3 IM0T,NTV»»E(2  1.31 
COMMON  /lO/  IDl.IDC 
DIMENSION  SUR SVE(42,3.4),REL( 2) 

this  routine  CALCULATES  RELAXED  STRESSES 


I 


1 


4V 


50 

52 

55 


60 


NSAVE»NTIMES( MISS) 

NEL*NT1M£S(M!  SS) 

00  50  1-2 .NUH 
00  50  J*2.3 

REL( 1 )«( ASUR/eSUR)*(CSUR**NFL-1.0) 

R6Li2)»<  ASuB/BSuB)»(CSU8**NFL-i.O) 

J1*J-I 

SIG«SURSVE  ( I . MISS.  J KREK  J1 ) 

CALL  INTERY(CUR^M.NMUN0.SURSVE( 1 .Miss ,1) , STR6S  S) 
Call  VLATEICURVM.SURSVE  II  .MISS.l ) .SIG.IVI  CL.NMCi^JO) 
NFL»NFl-2 

1F(NFl.E0«0>  GO  TO  A9 

ifiiviol.lt.o  ) go  to  I 

IF(S1G.lT.S1GL(3W>  go  TO  1 
NFL«NFL^1 

1FINFl.LT.NSA  yt)  NSAVe*K'FL 
IFINSAVE.EO.O  1 GO  TO  52 

continue 

NFl»N5AVE 

lFfMlSS.cO.NUMMlS»NTiMESCMISS)=NFL 
R£L1»I ASUR/BSOK)*(CSUR**NFL-l .» 

REl2«I ASU8/BSU6)*(CSUB**NFL-1 .» 

00  60  I-l.NUM 

SURSYEil.MlSS  .2)  = SURSVE(I.MISS,21-»RELl 

SURSVEl  l.MISS  .3)-SURSVEII .MISS.3»^REL2 

CO^^TINUE 

RETURN 

ENO 


1 

'i 
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